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Summary 


The cloud cover in a set of summertime and wintertime AVHRR data from the Arctic and 
Antarctic regions was analyzed using a pattern recognition algorithm. The data were collected 
by the NOAA-7 satellite on 6-13 January and 1-7 July 1984 between 60° and 90° north and 
south latitude in 5 spectral channels, at the Global Area Coverage (GAC) resolution of 
approximately 4 km. This data embodied a Polar Cloud Pilot Data Set which was analyzed by 
a number of research groups as part of a polar cloud algorithm intercomparison study. This 
study was intended to determine whether the additional information contained in the AVHRR 
channels (beyond the standard visible and infrared bands on geostationary satellites) could be 
effectively utilized in cloud algorithms to resolve some of the cloud detection problems caused 
by low visible and thermal contrasts in the polar regions. 

The analysis described here makes use of a pattern recognition algorithm which estimates 
the surface and cloud classification, cloud fraction, and surface and cloudy visible (channel 1 ) 
albedo and infrared (channel 4) brightness temperatures on a 2.5*x2.5* latitude-longitude grid. 
In each grid box several spectral and textural features were computed from the calibrated pixel 
values in the multispectral imagery, then used to classify the region into one of eighteen surface 
and/or cloud types using the maximum likelihood decision rale. A slightly different version of 
the algorithm was used for each season and hemisphere because of differences in categories 
and because of the lack of visible imagery during winter. The details of the classification 
algorithm are given in Ebert (1987a, 1987b, 1990). 

The classification of the scene is used to specify the optimal AVHRR channel for separating 
clear and cloudy pixels using a hybrid histogram-spatial coherence method. This method 
estimates values for cloud fraction, clear and cloudy albedos and brightness temperatures in 
each grid box. The choice of a class-dependent AVHRR channel allows for better separation 
of clear and cloudy pixels than does a global choice of a visible and/or infrared threshold. The 
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classification also prevents erroneous estimates of large fractional cloudiness in areas of cloud- 
free snow and sea ice. The hybrid histogram-spatial coherence technique and the advantages of 
first classifying a scene in the polar regions are detailed in Ebert (1989), which is also 
reproduced in the appendix of this report. 

The complete Polar Cloud Pilot Data Set has been analyzed and the results are presented 
and discussed in the next part of this report. The graphs were prepared with the assistance of a 
SSEC programmer, Mr. Frank LaFontaine. In April 1989 a magnetic tape containing the 
analyzed surface and cloud type, cloud fraction, clear and cloudy visible albedo and infrared 
brightness temperature in each 2.5°x2.5° box in every image, along with daily and weekly 
summary statistics, was forwarded to Professor E. Raschke at the University of Cologne for 
intercomparison with the analyses of other groups. The results may be obtained from Prof. 
Raschke at his new address at the GKSS, Hamburg. The results from the Arctic winter and 
Antarctic summer and winter analyses were presented at the IAMAP Symposium for Remote 
Sensing in the Polar Regions in Reading, England, in August 1989. A paper based on this 
presentation has been accepted for publication in the International Journal of Re mote Sensing, 
and the final manuscript is reprinted in the appendix. 
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Analysis of the Polar Cloud Pilot Data Set 

The data set contains 365 AVHRR images, or an average of 13 images per day. Each image 
was analyzed individually, then the analyses were combined into daily mosaics of classification 
and cloud fraction. These will be described in the following section for each season and 
hemisphere. The next section presents the weekly statistics of total, low and high cloud fraction. 
The fractional area and frequency of each surface and cloud type are given in the third section and 
compared to climatological values where possible. The last section examines the values of the 
surface and cloudy albedos and brightness temperatures for each of the surface and cloud classes. 

Reference maps of the Arctic and Antarctic are given in Figs. 1 and 2. 

1. Daily analyses 

This section describes in detail the nature and evolution of the cloud cover over each 7-day 
period. Daily mosaics of surface and cloud type were constructed by simply choosing, in each 
grid box, the classification which occurred most frequently during the day. The cloud fraction 
mosaics depict the daily mean total cloud fraction for the day. Thus it is possible that a region 
which has a clear classification in the mosaic may have a mean cloud fraction which is greater than 
zero if some cloud cover has been analyzed there during the day. Smoothing in the contouring 
routine may also produce this effect. 

Each mosaic has a "hole" immediately surrounding the pole where the imagery was not 
analyzed because of excessive distortion due to the earth's curvature. Also, some of the region 
between 90° E and 90°W (including the dateline) was not analyzed because the Pilot Data Set only 
included images which contained data in the hemisphere between 90*W and 90’E (including the 
Greenwich meridian). 

References are made to the synoptic weather during each period. (Charts for these days are 
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not shown here, but may be found in the publications of ECMWF, 1984, and the Free University 
of Berlin, 1984.) The goal is to relate the analyzed cloud cover to that which would be expected 
for various synoptic situations. 

a. Arctic summer 

Figures 3-9 show the automated classification and cloud fraction analyses for the Arctic 
during 1-7 July 1984 (days 183-189). The key to the surface and cloud types in the Arctic 
summer algorithm is given below: 


1 . Land without snow 

2. Land with snow 

3. Open water 

4. Unbroken sea ice 

5. Broken sea ice 

6. Ice sheet 

7. Mixed surface types 

8. Cs/As over land/water 

9. Ci/Cc/Ac over land/water 


10. St over land/water 

11. Sc over land/water 

12. Cu over land/water 

13. Cs/As over snow/ice 

14. Ci/Cc/Ac over snow/ice 

15. St over snow Ace 

16. Sc over snow/ice 

17. Cu over snow/ice 

18. High cloud over low cloud 


The standard meteorological abbreviations for cloud types are used (WMO, 1956). We assume 
that the low cloud types (classes 10, 11, 12, 15, 16, 17) consist of liquid water droplets, while 
the high cloud types (calsses 8, 9, 13, 14) are composed of ice crystals. Thus a mid-level cloud 
would be classified as either low or high, depending on its phase. Class 17, "cumulus over 
snow/ice,” includes isolated low clouds of small horizontal extent (i.e., much smaller than the grid 
square) since convective clouds over snow and ice are not common. "High clouds over low 
clouds” will frequently be abbreviated in the following sections as "layered clouds." The 
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classification skill for the Arctic summer training data was 84.6% (Ebert 1987b), and the root- 
mean-square (RMS) error in the cloud fraction was 0.18 when verified against a manual 
nephanalysis of the training imagery (Ebert 1989). 

Starting in the western part of the analyses, northwestern Canada showed significant cloud 
cover in the early part of the period associated with a decaying cyclone. These were replaced by 
high clouds in the latter part of the period. The Beaufort Sea had large amounts of low and 
layered clouds throughout most of the period, with clear patches of unbroken and broken sea ice 
on days 183 and 185-187. In the Canadian Archipelago layered clouds predominated, although 
there was a tendency for cirrus clouds to be found over Baffin Island. Baffin Bay had moderate 
amounts of stratus cloud until day 186 when high pressure moved in and significant clearing took 
place to reveal the broken sea ice cover. Moderate amounts of low, high, and layered clouds were 
analyzed over Greenland during the entire week. Clearer areas were found along the central 
eastern and western coasts on days 187-188. Stratus clouds covered the Greenland and 
Norwegian Seas, with fractional amounts generally exceeding 80%. Layered clouds on days 
188-189 were related to a trough connecting low pressure centers south of Iceland and northwest 
of Novaya Zemlya. This also produced large amounts of layered cloudiness in the Barents Sea 
and over Scandinavia, which was under low pressure during the entire week. Mainly clear 
conditions characterized Svalbard during the first two days of the period. Novaya Zemlya and the 
Kara Sea displayed moderate amounts of stratus and layered cloud cover. The analysis of "stratus 
over land/water" (class 10) on days 183-185 is erroneous, as the Kara Sea was neatly 100% ice- 
covered at that time. Siberia, particularly the eastern sector, remained mostly cloud-free during 
the week, with only small amounts of cumulus, stratocumulus, and cirrus clouds. 

b. Arctic winter 

Mosaics of classification and total cloud fraction for the Arctic during 6-12 January 1984 are 
shown in Figs. 10-16. The classification key for Arctic winter is: 
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1 . 

Land without snow 

10. St over land/water 

2. 

Land with snow 

1 1 . Sc over land/water 

3. 

Open water and thin ice 

12. Cu over land/water 

4. 

Unbroken sea ice 

13. Cs/As over snow/ice 

5. 

Broken sea ice 

14. Ci/Cc/Ac over snow/ice 

6. 

Ice sheet 

15. St over ice sheet 

7. 

Mixed surface types 

16. Sc over snow/ice 

8. 

Cs/As over land/water 

17. Inversion St over snow/ice 

9. 

Ci/Cc/Ac over land/water 

18. High cloud over low cloud 

Class 15, "stratus over ice sheet," designates a flat, fairly featureless cloud which occurred over 


Greenland. Class 17, "inversion stratus over snow/ice," refers to stratus clouds which are 
wanner than the underlying snow and ice surface. The Arctic winter algorithm had a classification 
skill of 71.9% and a RMS cloud fraction error of 0.24 on the training data (Ebert 1987a, 1990). 

Again starting in the west, high pressure over much of the Beaufort Sea during the week 
suppressed cloud cover, allowing only small amounts of stratus, stratocumulus, and layered 
clouds. On days 9-11 much of the clear area was classified as "broken sea ice" (class 5), when 
actually the sea ice cover was unbroken (class 4); the misclassification was probably the result of 
warm air advection. This illustrates one of the problems caused by the lack of visible data 
combined with relatively low spatial resolution. Stratus and layered clouds covered most of 
northwestern Canada until the last day of the period, while the Canadian Archipelago experienced 
primarily clear conditions. The frequent analysis of "ice sheet" (class 6) over Baffin and 
Ellesmere Islands reflects the cold textureless wintertime snow cover and should properly be 
labeled as "land with snow" (class 2) (although Ellesmere Island does have permanent ice cover at 
its higher elevations). Baffin Bay and Davis Strait also remained primarily clear throughout the 
period. Northern Greenland was cloud free through day 9, after which a large shield of high 
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clouds developed in the easterly flow and propagated slowly eastward. 

The conditions in southern Greenland, northern Europe, and the adjacent Greenland, 
Norwegian, and Barents Seas were dominated by the development and propagation of three 
cyclones. The first brought significant cloud cover with rain and snow to Sweden and Finland on 
days 7 and 8, while cumulus clouds filled the Norwegian Sea in the cold air behind the front. At 
this time a second cyclone was developing to the southeast of Greenland. As is intensified during 
days 8-10 the Greenland and Norwegian Seas were filled with layered and high clouds. The line 
of maximum cloud fraction along the occluded front can be seen in Figs. 14b and 15b. By day 12 
the low pressure center was decaying over Svalbard and the cloud cover was diminishing in the 
Norwegian and Barents Seas. The third cyclone moved into western Norway on day 12, bringing 
with it large amounts of layered and high clouds. 

The northwestern Soviet Union contained moderate amounts of cumulus and cirrus cloud 
cover early in the period, shifting to primarily stratiform cloudiness later on. Stratiform and high 
clouds over Novaya Zemlya and the Kara Sea on days 6 and 7 gave way to clear conditions 
during the next three days. Southwesterly flow ahead of the second cyclone led to significant 
high cloudiness on days 1 1 and 12. As in the summertime analysis, the stratiform clouds over the 
Kara Sea were misclassified as "stratus over land/water" (class 10) when in fact the surface was 
ice-covered. Sparse stratus and stratocumulus cloudiness was analyzed in the Laptev Sea, as well 
as in the central Arctic. However, the Arctic winter algorithm has difficulty detecting low clouds 
over sea ice at night (Table 3.12 of Ebert (1987a) indicates that "inversion stratus over snow/ice" 
was misclassified as clear 54% of the time) so the actual low cloud amounts in the central Arctic 
are probably much higher. In Siberia, cirrus and stratus clouds were present in large quantities 
between 60°E and lOO’E during the entire period, probably as the result of moist air advection 
from the west. Further to the east clearer conditions prevailed. 
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c. Antarctic summer 

Figures 17-23 show the cloud analyses for the Antarctic on 6-12 January 1984. The analyses 
are frequently incomplete near the northern edges because the satellite crossed 60*S in the early 
morning hours when part of the image was still dark. Since the Antarctic summer algorithm relies 
on data in the visible channels, the affected pixels were discarded. The classes used in the 
Antarctic summer algorithm are: 


1 . Ice haze 

2. Land with snow 

3. Open water 

4. Unbroken sea ice 

5. Broken sea ice 

6. Ice sheet 

7. Mixed surface types 

8. Cs/As over land/water 

9. Ci/Cc/Ac over land/water 


10. St over land/water 

11. Sc over land/water 

12. Cu over land/water 

13. Cs/As over snow/ice 

14. Ci/Cc/Ac over snow/ice 

15. St over snow/ice 

16. Sc over snow/ice 

17. Cu over snow/ice 

18. High cloud over low cloud 


Class 1 differs from the Arctic summer algorithm. A cloud type characteristic of low level ice 
crystal haze had a tenuous appearance in the visible imagery, a small signal in channel 3, and a 
slightly warmer brightness temperature than the surface in channel 4. The classification skill for 
the Antarctic summer algorithm was 85.1% when tested on the training data (Ebert 1987a, 1990). 

In general, Antarctica was characterized by quite small fractional cloud cover over the interior 
of the continent, moderate amounts near the coasts, and nearly overcast conditions in the 
surrounding ocean. Cyclones developing and propagating eastward in the Southern Ocean near 
60*S tended to have large amounts of low and layered clouds in their centers and eastern flanks. 
Examples of this are seen at 10'W on day 6 (see Fig. 3b of Ebert, 1990, for the sea level pressure 
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chart for this day), 30' W on day 10, and especially 1 10'E on day 11, where the central pressure 
of 957 mb and the strong upper level support indicated a mature cyclone. The Weddell Sea and 
Ronne Ice Shelf were persistently covered by extensive stratus and layered cloud cover. The 
region just to the west of the Antarctic Peninsula was less cloudy. This is unusual from a climatic 
point of view, but the easterly air flow over the peninsula between days 6 and 10 may have caused 
downslope subsidence warming which suppressed the cloud cover. 

Western Antarctica had varying amounts of low clouds during the week, with stratus as the 
dominant type. Layered clouds occurred close to the coast and near the base of the Antarctic 
Peninsula. Clouds occurred frequently over East Antarctica as well, but their horizontal extent 
was typically small so the total cloud fraction was often less than 20%. Low cloud types were the 
most common except in regions where the upper tropospheric circulation brought moister air from 
the north to form cirrus clouds (for example, at 45'E on day 8 and at 140'E on day 11). 

d. Antarctic winter 

The Antarctic classification and cloud fraction mosaics for 1-7 July 1984 (days 183-189) are 
shown in Figs. 24-30. Again, the analysis is incomplete near 60°S, this time because of solar 
radiation contaminating some of the channel 3 radiances. Since the wintertime algorithm requires 
purely thermal radiances, the affected pixels were discarded. The classification key for the 
Antarctic winter algorithm is given below: 

1. Polar stratospheric cloud 10. St over land/water 

2. Land with snow 11. Sc over land/water 

3. Open water and thin ice 12. Cu over land/water 

4. Unbroken sea ice 13. Cs/As over snow/ice 

5. Broken sea ice 14. Ci/Cc/Ac over snow/ice 

15. St over snow/ice 

16. Sc over snow/ice 


6. Ice sheet 

7. Mixed surface types 



8. Cs/As over land/water 

9. Ci/Cc/Ac over land/water 


17. Inversion St over snow/ice 

18. High cloud over low cloud 
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The Antarctic winter algorithm had a classification skill of 78.4% on the training data (Ebert 
1987a, 1990). 

The significant feature in this set of imagery was the large 2-lobed polar stratospheric cloud 
over East Antarctica during the entire week. The brightness temperatures in the center of the lobes 
reached as low as 179 K, a full 19 degrees colder than the surface temperature recorded at nearby 
Vostok station. The calibration error of the AVHRR channel 4 sensor probably does not exceed a 
few degrees (see Appendix A of Ebert, 1990), so one plausible explanation for the very cold 
temperatures is the presence of an optically thick cloud at high altitude. The boundaries of this 
cloud are not sharp in the imagery-the sharp boundaries seen in the cloud fraction analysis are an 
artifact of the algorithm's temperature differentiation between "polar stratospheric cloud" (class 1), 
"cirrostratus over snow/ice" (class 13), and "ice sheet" (class 6), which all look qualitatively 
similar. The distinction between areas of "polar stratospheric cloud" and "cirrostratus over 
snow/ice” is in fact rather arbitrary and perhaps need not be made. 

Qoud amounts were generally small over the rest of the continental interior, with inversion 
stratus being the dominant type. Even so, the true fractional coverage of inversion stratus may be 
greater because the Antarctic winter algorithm has difficulty detecting this cloud type over snow 
and ice at night (Table 3.15 of Ebert (1987a) indicates that 38% of class 17 training samples were 
misclassified as clear). Layered cloudiness occurred in greater amounts near the coastlines and in 
the vicinity of cyclones. Good examples of this are found near 120*W on days 184-188 and near 
90"E throughout the period, where a strong ridge at 120‘E blocked the movement of cyclones 
around the periphery of the continent. Stratus, stratocumulus, and layered clouds were analyzed 
in the Bellingshausen and Amundsen Seas and between 0* and 60*E. Cloud-free unbroken sea ice 
was analyzed in the northern Weddell Sea, while small amounts of stratus and occasional cirrus 
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clouds were found closer to the Ronne Ice Shelf. The Weddell Sea constituted an outflow region 
for cold, dry continental air. 

The automated analyses have not been independently verified against surface observations 
except for the first day of each analysis period (Ebert 1989, 1990). Such a validation encounters 
significant problems and uncertainties with regard to spatial and temporal resolution, not to 
mention the extreme sparsity of surface observations at the highest latitudes. However, the 
analyses described above show a high degree of coherence from day to day, as would be expected 
for naturally evolving and propagating systems. Also, the analyzed cloud types and amounts are 
consistent with the synoptic weather situations given by the ECMWF and Free University of 
Berlin charts. These two factors give credence to the algorithm and its results. 

2. Weekly analyses 

The mosaics of weekly mean total, low, and high cloud cover were averaged from all of the 
analyses during each week. The standard deviation of the total cloud fraction is a measure of the 
temporal variability, where low values correspond to persistent conditions. The mean low and 
high cloud fractions give a rough idea of the vertical distribution of cloud cover. The low cloud 
types include stratus, inversion stratus, stratocumulus, cumulus, ice haze, and the low cloud 
portion of class 18, "high cloud over low cloud." The high cloud types were 
cirrus/cirrocumulus/altocumulus (classes 9 and 14), cirrostratus/altostratus (classes 8 and 13), 
polar stratospheric cloud, and the high cloud portion of "high cloud over low cloud.” No 
corrections for cloud overlap have been made to account for low clouds hidden from the view of 
the satellite, so the low cloud fraction analyzed for class 18 is likely to significantly underestimate 
the true low cloud amount. 

For comparison, the high latitude climatologies of Huschke (1969) and Schwerdtfeger (1970) 
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are shown in Fig. 31. These are based on relatively few surface observations and reconnaissance 
flights and may not be extremely reliable. (There is some feeling in the polar community that the 
Huschke summertime cloud amounts may be too high (Curry and Herman 1985, McGuffie et al. 
1988).) However, they do convey a general sense of the annual cycle of cloudiness in the Arctic 
and Antarctic. The cloud fraction in the Arctic is much smaller in winter than in summer, and the 
Canadian Arctic tends to be less cloudy than the other regions. The Antarctic has a less 
pronounced annual cycle, but a greater contrast between the extensive cloudiness over the ocean 
and clearer conditions inland. 

a. Arctic summer 

Figures 32 a and b show the mean and standard deviation of total cloud fraction for the week 
of 1-7 July 1984. Figures 32 c and d separate the total cloud fraction into its low and high cloud 
components. The clearest regions were over land in Siberia and in the Canadian Archipelago, 
while the maximum cloud cover during the week was found over the Norwegian and Barents 
Seas. The greatest persistence was found here, with standard deviations below 20%. The 
greatest variability occurred in northern Canada in the vicinity of the decaying cyclone, and in the 
eastern Arctic. Low clouds occurred much more frequently than high clouds everywhere except 
in northwestern Canada where cirrus and layered clouds predominated. Large amounts of low 
cloudiness are typical of the summertime Arctic (Vowinckel and Orvig 1970). However, the 
analyzed total cloud fraction was everywhere smaller than the Huschke (1969) climatological 
values for July. 

b. Arctic winter 

The weekly mosaics for the Arctic during 6-12 January 1984 are shown in Fig. 33. Overall, 
the total cloud cover was less than in the summer case, being below 20% over much of the 
Canadian and Siberian Arctic. These amounts are far less than the climatological amounts of 
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Huschke (1969), partly due to an underestimation of stratus amounts as discussed in Section lb. 
The cloudiest regions occurred over northwestern Canada, which had persistent layered clouds, 
and in northwestern Siberia. Vigorous synoptic activity led to high overall variability in the total 
cloud fraction, with standard deviations exceeding 40% in much of the Arctic. In contrast to the 
summer case, analyzed high cloud amounts were greater than low cloud amounts everywhere but 
in the Siberian Arctic. However, this should be interpreted with caution because the satellite 
cannot detect low clouds beneath a higher cloud layer. Also, many mid-level clouds of ice phase 
were grouped with high clouds in winter, while summertime mid-level clouds were more likely to 
be composed of liquid water and grouped with low clouds. 

c. Antarctic summer 

Figure 34 shows the weekly mosaics for the Antarctic on 6-12 January 1984. The total cloud 
fraction resembled the daily mosaics in that the continental interior was mostly cloud-ffee while 
the Southern Ocean had large fractional cloud amounts. This resembles the climatological field in 
January (Fig. 31b). The minimum standard deviation was less than 10% in the Weddell Sea, 
which was characterized by persistent stratus and layered clouds. The greatest variability occurred 
near the coasts of the continent where synoptic cyclones had a strong influence on the cloud 
conditions. Low cloud fractions were greater than high cloud fractions in all regions of the 
Antarctic except along the coast at 150°E. This is partially a result of warmer summertime 
temperatures favoring water phase clouds, as discussed above. 

d. Antarctic winter 

The weekly Antarctic mosaics for the week of 1-7 July 1984 are given in Fig. 35. The two- 
lobed polar stratospheric cloud and the clear regions on the Antarctic Plateau and in the Weddell 
Sea are markedly evident on the mean total cloud fraction map. The variabilities were also quite 
low in these regions, with standard deviations ranging between 10% and 20%. Large standard 
deviations were analyzed at the edges of the polar stratospheric cloud where slight shifts in 
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position and temperature caused a shift in the analyzed sharp cloud boundary (refer to Section Id). 
Except for the Weddell Sea, the oceanic regions showed moderate mean cloudiness with a 
standard deviation of approximately 30%, reflecting the evolution and motion of clouds during the 
7-day period. This is somewhat smaller than the climatological cloud fraction in July given by 
Schwerdtfeger (1970) and may again be related to difficulties in analyzing low clouds over sea ice 
at night. Bearing in mind this shortcoming, low clouds and high clouds were present in roughly 
equal proportions in the Bellingshausen and Amundsen Seas and between 0” and 60'E. Low 
clouds were more common in the southern Weddell Sea, while high cloud predominated at the tip 
of the Antarctic Peninsula, in the ocean east of 60°E, and of course in the locations of the polar 
stratospheric cloud. 

3. Surface and cloud type statistics 

Figures 36-39 show the fractional area statistics of two surface types and five cloud types for 
each analysis period. In this analysis the layered clouds (class 18) have been assumed to contain 
Ci/Cc/Ac in the upper layer and stratus in the lower layer. The fractional area covered by each 
surface or cloud type is a quantitative measure of the prevalence of that type. The trends in the 
curves are a result of synoptic evolution during the week. 

The mean frequency indicates the fraction of grid boxes which were classified as each surface 
and cloud type, and indicates how often that type occurred during the given period. Table 1 gives 
values of frequency for the four analysis periods. One can compare this to similar statistics 
computed by Hahn et al. (1982, 1984). Their cloud frequency represents the fraction of all ship 
and ground observations reporting each cloud type over a span of twelve years. Their values for 
ocean and land were area weighted by give zonal mean values; these are summarized in Table 2, 
and will be referred to as the "observed climatological frequencies” in the following discussion. 
The comparison is inexact because the two quantities differ with respect to cloud views (from 
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space versus from the surface), spatial extent (60 -90 N versus 60 -80 N, also 2.5 x2.5 grid 
boxes versus point observations *), temporal domain (specific 7-day periods versus seasonal 
means), and cloud categories (compare Tables 1 and 2). Nevertheless, such a comparison is 
instructive because it indicates departures of the analyzed polar cloud fields of 6-12 January and 
1-7 July from climatology. 

a. Arctic summer 

The mean cloud fraction in the Arctic for the week of 1-7 July 1984 was 0.56, with only small 
daily departures from this value. This is in contrast to Huschke s (1969) climatological value of 
0.80 in July (however, his analysis does not include land regions). The area of clear land and 
water was much greater than that for clear snow and ice, first because the sea ice extent was 
nearing its summertime minimum, and also because the cloud-free conditions in Siberia had a 
large influence on the mean value. The frequency of grid boxes analyzed as cloud-free was 0. 15, 
compared to the observed climatological value of 0.05 between 60* and 80*N. Stratus was by far 
the most extensive cloud type analyzed, and was diagnosed in 61% of the grid boxes. This was 
very close to the observed climatological frequency of 0.60 for combined stratus, stratocumulus, 
and nimbostratus. The fractional area covered by cumulus clouds was less than 1%. High 
stratiform clouds were less common than broken high clouds. Their combined frequency of 0.39 
was on the same order as the observed climatological frequency of 0.45. 

b. Arctic winter 

The analyzed total cloud fraction of 0.37 in the Arctic was much smaller in the winter period 
than in the summer period, and smaller than the Huschke (1969) January climatological cloud 
fraction of 0.50. The large amount of clouds created by the passage of cyclones over Europe and 
Western Siberia was offset by the large clear regions in the central Arctic and Canadian 


1 The Hahn et al. (1982, 1984) statistics are also biased by the locations of the surface stations and shipping 
lanes. 
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Archipelago. The area of clear snow and ice was much greater than in summer, due partly to the 
seasonal increase in sea ice extent, but more importantly to the significant reduction in stratus 
cloud cover. As the total amount of high cloudiness increased in the Barents and Chukchi Seas, 
the area of clear snow and ice decreased (Fig. 37). The mean frequency of high clouds (0.21) 
was less than the mean frequency of low clouds (0.28), but their mean fractional coverage was 
greater (0.23 as opposed to 0.15), at least when viewed from above. The analyzed high cloud 
frequency was significantly lower than the observed climatological frequency of Ci/Cs/Cc/As/Ac 
of 0.56. The analyzed frequency for stratus and stratocumulus of 0.23 was also much lower than 
the observed stratiform frequency of 0.47, but some low clouds were undoubtedly missed by the 
classification algorithm. The week of 6-12 January 1984 may have been less cloudy than occurs 
climatologically. 

c. Antarctic summer 

The mean cloud fraction analyzed in the Antarctic region for the week of 6-12 January 1984 
was 0.45, a balance between the very cloudy conditions in the Southern Ocean and the sparse 
cloudiness in the interior. The cloud cover was composed primarily of stratus and citrus clouds, 
with cirrostratus and cumulus clouds occurring only rarely. There were only four grid boxes 
analyzed as clear ocean, for a frequency of 0.0001. The observed climatological frequency of 
clear sky is 0.01 between 60°S and 70°S. The analyzed frequency of 0.60 for combined stratus 
and stratocumulus cloudiness was comparable to the observed climatological frequency of 0.64 
for combined St/Sc and Ns clouds in the Antarctic. However, the computed cirrus frequency of 
0.18 was much smaller than the observed climatological value of 0.48. This may partially explain 
why the total cloud fraction was lower than has been climatologically observed (refer to Fig. 31b). 

d. Antarctic winter 

The Antarctic mean cloud fraction of 0.38 during 1-7 July 1984 comprised roughly equal 
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amounts of cirrus, stratus, and polar stratospheric cloud, with much smaller quantities of 
cirrostratus, stratocumulus, and cumulus. The true amount of stratus cover may have been 
underestimated. In spite of that, stratus was the most commonly analyzed cloud type, with a 
frequency of 0.43. The data used by Hahn et al. (1982) did not contain enough cloud 
observations over the Southern Ocean to derive reliable cloud frequencies. They did, however, 
obtain enough clear observations to derive an observed climatological clear-sky frequency of 0.25 
between 6(TS and 90°S. The analyzed clear-sky frequency of 0.23 compares closely to this 
value. The cumulus frequency of 0.001 was the lowest of the cloud types. 

4. Surface and cloudy albedos and brightness temperatures 

In this section we examine the values of channel 1 (0.6 |im) albedo and channel 4(11 pm) 
brightness temperature for the clear and cloudy portions of each of the eighteen surface and cloud 
classes (denoted a g T s , a c , and T c respectively). These were derived using the hybrid 
histogram-spatial coherence method. 

The measured visible reflectances were normalized by the cosine of the solar zenith angle, but 
no other corrections for atmospheric effects were made. The term "cloudy" is used instead of 
"cloud" because the effects of the underlying surface on the measured albedo and brightness 
temperature, namely multiple reflections in channel 1 and upwelling longwave radiation in channel 
4, have not been removed. Thus a semi-transparent cloud appears brighter and colder over a 
snow surface than over a land surface. This is seen in the results to be presented here, although 
there may also be morphological differences between clouds over land and clouds over snow. In 
principle it is possible to derive cloud optical properties from the satellite measured albedos and 
brightness temperatures given appropriate microphysical models (Kriebel et al. 1989), and to use 
these to estimate cloud temperature and surface radiative fluxes. (In fact, that was one of the 
original reasons for classifying cloud types when this study was begun.) In practice, however, 
small errors in analyzed cloud fraction and surface and cloudy albedos and brightness 
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temperatures, combined with deviations of the true cloud microphysical properties from the 
idealized models, led to large errors in the cloud optical properties and unreasonable values and 
spatial variations in cloud temperatures and surface fluxes. These difficulties are certainly greater 
in the polar regions than elsewhere because a high proportion of clouds are semi-transparent and 
because multiple reflections are prevalent. 

The means and standard deviations of surface and cloudy albedos and brightness temperatures 
are presented in both tabular and graphical form for each of the four analysis periods. 

a. Arctic summer 

Table 3 and Fig. 40 show the mean and standard deviation of a s T s , a c , and T c for 1-7 July 
1984. Open water had the darkest surface, with a s =0.05, and clear land had an albedo of 
a s =0.09. The albedo of unbroken sea ice, much of which was in a state of melting, was 0.49, 
equal to the 0.6 |Xm spectral albedo of an Arctic early season melt pond measured by Grenfell and 
Maykut (1977). Clear snow was the brightest surface in the imagery, with a s between 0.61 and 
0.65 (classes 2, 6, 17). The mean surface albedo for class 18 had an intermediate value of 0.34 
and a large standard deviation because layered clouds were analyzed over many different surface 
types. The albedos for clouds occurring over snow and ice were much higher than those for 
clouds over land and water. The cumulus cloudiness over snow and ice had the highest mean 
albedo of 0.70. 

Channel 1 separated surface and clouds quite well in classes 1-12. The channel 4 brightness 
temperature was effective for separating classes 8, 9, 11, 12, 13, and 18. (Classes 15-18 used 
channel 3 as the class-optimal channel but channel 3 radiances are not shown here.) The land had 
the warmest temperatures, at T s =235-290 K (classes 1, 2, 12). The brightness temperature of 
water was approximately 277 K (classes 3, 5). The mean measured brightness temperature of sea 
ice was 271.3 K, close to the melting point of sea ice of 271.2 K (Untersteiner 1961). Clouds 
over snow and ice tended to have colder brightness temperatures and smaller variabilities than 
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those over land and water. The exception was "Cs/As over land/water," which had a mean T c of 
243.8 K, but with a large variability. The brightness temperatures of the more tenuous cirrus 
clouds were 10-20* greater than that of cirrostratus. Cumulus clouds over land and water had a 
mean value of T c =272.6 K. These cloud appeared to have small horizontal and vertical 
dimensions and were most frequently found over Siberia. 

b. Arctic winter 

The means and standard deviations of surface and cloudy channel 4 brightness temperatures 
are given in Table 4 and Fig. 41. Because of the high level of noise in channel 3 and the near 
redundancy of channel 5 with channel 4, the class-optimal channel was chosen as channel 4 for all 
classes. 

All occurrences of cloud-free land were diagnosed as being snow-covered with a mean 
brightness temperature of 270.9 K. Open water and thin ice was only slightly warmer at 
T s =272.3 K. The surface brightness temperature analyzed for classes 8-12 had high standard 
deviations because they included both the relatively warm Norwegian Sea and the cold land of 
Canada and Siberia (cloud over land/water included cloud overlying snow-covered land in this 
analysis). Similarly the surface temperatures for classes 13-17 had high variabilities because they 
represent both sea ice (T s =233-247 K) and ice sheet (T s =225 K). The unbroken sea ice statistics 
include some misclassified cases of inversion stratus, and may be slightly too warm. The low 
standard deviations of T s and T c computed for class 17 may have resulted from the algorithm 
correctly analyzing the "obvious" cases while failing to recognize the less obvious stratus cases. 
The differences between surface and cloudy brightness temperatures in classes 8-12 were 
generally sufficient to enable them to be distinguished using channel 4. This was less true for 
classes 13-18, and was undoubtedly a large source of error in estimating the fractional cloudiness 
over snow and ice in the wintertime Arctic. 
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c. Antarctic summer 

Means and standard deviations of a s T s , a c , and T c for the Antarctic during 6-12 January 
1984 are given in Table 5 and Fig. 42. Most of the comments made for the Arctic summer case 
apply qualitatively to the Antarctic summer results. The snow albedo was slightly brighter in the 
Antarctic (a s =0.62-0.70), as was the mean sea ice albedo (a s =0.60) because surface melting is 
much less pronounced in the southern polar region. The mean surface albedo of ct s =0. 1 0 for 
class 18 is certainly erroneous because layered clouds were analyzed over all surface types. ^ 
Channel 1 effectively separated clouds from land and water, but was quite useless for 
distinguishing clouds from snow and ice over Antarctica. In fact, some of the cloudy albedos 
were actually lower than the snow albedos. This can happen when the cloud itself is less 
reflective than the surface, and also if some of the shadowed pixels are included in the analysis as 
was the case for ice haze. 

The water temperature in the Antarctic (T s =272 K) was much colder than in the Arctic were 
warm currents influence the temperature of the Norwegian Sea. The brightness temperature of 
snow and ice ranged from 240.4 K (class 1) on the high plateau to 268.3 K for unbroken sea ice. 
The surface elevation strongly influences the snow and cloud temperatures. This is why 
brightness temperatures for low cloud classes (15, 16, 17), which were diagnosed over all parts 
of the continent, were colder than those for the high cloud classes (13, 14), which were diagnosed 
close to the coast. The standard deviation for T s for classes 13-17 was also high. T s for class 18 
was too high for the reasons explained in Section 4a, and probably led to the overestimation of 
layered cloud fraction in West Antarctica. The high cloud brightness temperatures over land a 
water (classes 8, 9) were slightly warmer than the low cloud brightness temperatures because they 
were much more tenuous. 


2 After first analyzing the grid boxes closest to 60°S and building a class characteristic value of a s appropriate for 
the ocean surface, when the algorithm encountered layered clouds over snow and ice it "forced" a s to be close to the 
class characteristic value. The specific details of this method are described in Ebert (1989). This obvious defect is 
being corrected in an updated version of the algorithm. 
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d. Antarctic winter 

Table 6 and Fig. 43 show the means and standard deviations of T s and T c for the week of 1-7 
July 1984 in the Antarctic. Because the ocean surface was covered with sea ice in most of the 
domain, there were few incidences of open water and this ice (class 3) or clouds over land/water 
(classes 8-12). The mean brightness temperature of open water and thin ice in cloud-free 
conditions was 267.6 K, while beneath clouds it tended to be warmer by a few degrees. This is 
physically reasonable if the radiative loss at the surface under clear skies led to rapid cooling and 
refreezing. As in the Antarctic summer, the surface temperature beneath clouds of classes 13-17 
had a large range and variability due to the effects of surface elevation. The brightness 
temperature of the ice sheet ranged between 212 K and 229 K (classes 1, 6, 13). The cloudy 
brightness temperatures for clouds over land and water were 20-35 greater than those for 
corresponding cloud types over snow and ice. The coldest clouds were polar stratospheric cloud 
(T c =199 K) and cirrostratus over snow and ice (T c =203 K). As discussed earlier, the distinction 
between these two cloud types may not be meaningful. Using channel 4 alone it was possible to 
separate cloud from surface in most classes except 15 and 17, which designate stratus over snow 


and ice. 
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Table 1. Frequencies of surface and cloud types. 
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Table 2. Seasonally and zonally averaged frequencies of cloud types based on data from Hahn et al. (1982, 1984). The quantity 
outside the parentheses is the area weighted frequency; the quantites inside the parentheses are the frequencies over ocean and over 
land. 



Surface Cloudy Surface Ch. 4 Cloudy Ch. 4 

Number Class- Ch. 1 Ch. 1 Brightness Brightness 

ofObser- optimal Albedo Albedo Temperature (K) Temperature (K) 

Class vations Channel a. a T s T s 
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Surface Ch. 4 Cloudy Ch. 4 
Class- Brightness Brightness 

Number of optimal Temperature (K) Temperature (K) 
Class Observations Channel T s T c 


1. Land 0 

2. Snow on land 1119 

3 . Open water and thin ice 267 

4. Unbroken sea ice 9202 

5 . Broken sea ice 1 2557 

6. Ice sheet 1529 

7. Mixed surface types 2 2230 

8 . Cs/As over land/water 298 

9. Ci/Cc/Ac over land/water 2897 

10. St over land/water 2010 

11. Sc over land/water 272 

12. Cu over land/water 1680 

13. Cs/As over snow/ice 1622 

14. Ci/Cc/Ac over snow Ace 374 

15. St over ice sheet 551 

16. Sc over snow/ice 2937 


1 7 . Inversion St over snow/ice 3003 

18. High cloud over low cloud 2390 


4 


4 

270.9 

(3.5) 



4 

272.3 

(2.4) 



4 

233.7 

(7.0) 



4 

270.0 

(0.4) 

243.2 

(4.2) 

4 

225.0 

(5.6) 



4 

260.9 

(5.1) 

237.2 

(7.6) 

4 

261.3 

(3.2) 

227.3 

(6.1) 

4 

261.1 

(4.7) 

245.4 

(5.5) 

4 

260.4 

(3.7) 

246.9 

(6.9) 

4 

265.1 

(7.5) 

243.6 

(4.6) 

4 

267.4 

(7.1) 

250.2 

(3.1) 

4 

232.7 

(6.9) 

227.0 

(4.7) 

4 

246.5 

(9.9) 

231.4 

(8.7) 

4 

236.9 

(6.3) 

221.3 

(5.9) 

4 

241.3 

'4.6) 

231.3 

(3.6) 

4 

232.4 

U-3) 

242.5 

(3.7) 

4 

248.6 

(9.7) 

235.0 

(7.4) 


Table 4. Means (standard deviations) of surface and cloudy brightness temperatures for the 
Arctic during 6-12 January 1984. 


lH Surface" refers to open water, "Cloud” refers to sea ice. 
2 "Surface" refers to warmer surface, "Cloud" to colder surface. 


Surface Cloudy Surface Ch. 4 Cloudy Ch. 4 

Number Class- Ch. 1 Ch. 1 Brightness Brightness 

of Obser- optimal Albedo Albedo Temperature (K) Temperature (K) 

Class vations Channel a* ' a„ T c T c 
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^’Surface" refeis to land, "Cloud" refers to snow. 

2 "Surface" refers to open water, "Cloud" refers to sea ice. 

3 "Surface" refers to warmer, daiker surface, "Cloud" to colder, brighter surface. 
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Surface Ch. 4 Cloudy Ch. 4 
Class- Brightness Brightness 

Number of optimal Temperature (K) Temperature (K) 
Class Observations Channel T s T c 


1 . Polar stratospheric cloud 6687 

2. Snow on land 18 

3 . Open water and thin ice 74 

4. Unbroken sea ice 2444 

5 . Broken sea ice 1 53 

6. Ice sheet 5698 

7 . Mixed surface types 2 267 

8 . Cs/As over land/water 3 

9. Ci/Cc/Ac over land/water 41 

10. St over land/water 56 

11. Sc over land/water 606 

12. Cu over land/water 55 

13. Cs/As over snow/ice 1456 

14. Ci/Cc/Ac over snow/ice 2898 

15. St over snow/ice 6191 

16. Sc over snow/ice 795 


17. Inversion St over snow/ice 6906 

1 8 . High cloud over low cloud 3052 


4 

211.8 

(3.6) 

199.0 

(6.4) 

4 

244.5 

(5.3) 



4 

267.6 

(1.9) 



4 

242.7 

(5.9) 



4 

268.0 

(2.1) 

249.2 

(3.2) 

4 

220.0 

(8.2) 



4 

258.1 

(7.7) 

238.7 

(6.4) 

4 

267.1 

(0.6) 

237.8 

(2.9) 

4 

270.0 

(2.1) 

244.6 

(5.4) 

4 

268.1 

(2.0) 

253.5 

(4.3) 

4 

269.8 

(1.9) 

260.5 

(4.2) 

4 

271.5 

(1.4) 

251.8 

(6.3) 

4 

228.8 

(10.8) 

202.8 

(7.8) 

4 

241.9 

(7.6) 

224.9 

(4.6) 

4 

246.8 

(4.9) 

234.9 

(5.7) 

4 

254.4 

(7.5) 

238.8 

(4.2) 

4 

223.4 

(9.4) 

244.6 

(4.2) 

4 

266.8 

(0.9) 

237.8 

(12.7) 


Table 6. Means (standard deviations) of surface and cloudy brightness temperatures for the 
Antarctic during 1-7 July 1984. 


^’Surface" refers to open water, "Cloud" refers to sea ice. 
2 ”Surface" refers to wanner surface, "Cloud" to colder surface. 


Figure Captions 
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Fig. 1. The north polar region. 

Fig. 2. The south polar region, including some manned research stations. 

Figs. 3-9. (a) Automated classification and (b) cloud fraction analyses for the Arctic during 
1-7 July 1984. Refer to p. 4 for a key to the surface and cloud types for the Arctic 
summer algorithm. 

Figs. 10-16. (a) Automated classification and (b) cloud fraction analyses for the Arctic 
during 6-12 January 1984. Refer to p. 5 for a key to the surface and cloud types for the 
Arctic winter algorithm. 

Figs. 17-23. (a) Automated classification and (b) cloud fraction analyses for the Antarctic 
during 6-12 January 1984. Refer to p. 7 for a key to the surface and cloud types for the 
Antarctic summer algorithm. 

Figs. 24-30. (a) Automated classification and (b) cloud fraction analyses for the Antarctic 
during 1-7 July 1984. Refer to p. 9 for a key to the surface and cloud types for the 
Antarctic winter algorithm. 

Fig. 31. Annual variation of cloud fraction with region in (a) the Arctic (Huschke 1969), 
and (b) the Antarctic (from the climatic tables of Schwerdtfeger 1970). 

Fig. 32. (a) Mean total cloud fraction, (b) standard deviation of total cloud fraction, (c) 
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mean low cloud fraction, and (d) mean high cloud fraction for the Arctic during 1-7 July 
1984. 


Fig. 33. As in Fig. 32, for the Arctic during 6-12 January 1984. 

Fig. 34. As in Fig. 32, for the Antarctic during 6-12 January 1984. 

Fig. 35. As in Fig. 32, for the Antarctic during 1-7 January 1984. 

Fig. 36. Fractional area covered by each surface and cloud type during the week of 1-7 
July 1984 in the Arctic. 

Fig. 37. Fractional area covered by each surface and cloud type during the week of 6-12 
January 1984 in the Arctic. 

Fig. 38. Fractional area covered by each surface and cloud type during the week of 6-12 
January 1984 in the Antarctic. 

Fig. 39. Fractional area covered by each surface and cloud type during the week of 1-7 
July 1984 in the Antarctic. 

Fig. 40. Mean and standard deviation of channel 1 albedo and channel 4 brightness 
temperatures for the Arctic during 1-7 July 1984. Closed circles denote the surface 
values, open circles denote the cloud values, and error bars indicate one standard 
deviation on each side of the mean. For the interpretation of "cloud" in classes 2, 5, and 
7, refer to the footnote in Table 3. 
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Fig. 4 1 . Mean and standard deviation of channel 4 brightness temperatures for the Arctic 
during 6-12 January 1984. For the interpretation of "cloud" in classes 5 and 7, refer to 
the footnote in Table 4. 

Fig. 42. Mean and standard deviation of channel 1 albedo and channel 4 brightness 
temperatures for the Antarctic during 6-12 January 1984. For the interpretation of 
"cloud" in classes 2, 5, and 7, refer to the footnote in Table 5. 

Fig. 43. Mean and standard deviation of channel 4 brightness temperatures for the 
Antarctic during 1-7 July 1984. For the interpretation of cloud in classes 5 and 7, refer 


to the footnote in Table 6. 
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Figure 31 . Annual variation of cloud fraction with region in (.a) 
the Arctic (Huschke, 1969), and (b) the Antarctic (from the 
climatic tables of Schwerdtfeger, 1970). 
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Fig. 36. Fractional area covered by each surface and cloud type during the week of 1-7 July 1984 in the Arctic. 




Arctic 6-12 January 1984 
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Fig. 37. Fractional area covered by each surface and cloud type during the week of 6-1 2 January 1984 in the Arctic. 




Antarctic 6-12 January 1984 
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Fig. 38. Fractional area covered by each surface and cloud type during the week of 6-12 January 1984 in the 
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Fig. 39. Fractional area covered by each surface and cloud type during the week of 1-7 July 1984 in the Antarctic. 




Arctic summer 


76 




o 

o 

00 


o 

CM 


<D 


O 3 
00 “E 


a> 

Q. 

E 

CM 0) 


O 


O 

CO 

CM 


CO 

0) 

0 

c 


o 

CO 

CM 


O) 

CD 


o 

rj- 

CM 


00 

d 


CO 

d 


O 

Tf X3 

• CD 


CM 

d 


o 

o 


Fig. 40. Mean and standard deviation of channel 1 albedo and channel 4 brightness temperature for the Arctic during 1-7 July 1984. 
Closed circles denote the surface values, open circles denote the cloud values, and error bars indicate one standard deviation on 
each side of the mean. For the interpretation of "cloud" in classes 2, 5, and 7, refer to the footnote in Table 3. 
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Fig. 41 . Mean and standard deviation of channel 4 brightness temperature for the Arctic during 6-12 January 1984. 
For the interpretation of "cloud" in classes 5 and 7, refer to the footnote in Table 4. 
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Fig. 42. Mean and standard deviation of channel 1 albedo and channel 4 brightness temperature for the Antarctic during 6-13 
January 1984. For the interpretation of "cloud" in classes 2, 5, and 7, refer to the footnote in Table 5. 
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Fig. 43. Mean and standard deviation of channel 4 brightness temperature for the Antarctic during 1-7 July 1984. 
For the interpretation of "cloud" in classes 5 and 7, refer to the footnote in Table 6. 
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Two papers based on the analysis of the Polar Cloud Pilot Data Set are included in this 
appendix: 

Ebert, E.E., 1989: Analysis of polar clouds from satellite imagery using pattern 

recognition with a statistical cloud analysis scheme. J. Appl. Meteor., 28, 382-399. 

Ebert, E.E., 1990: Pattern recognition analysis of polar cloud during summer and winter. 
Int. J. Remote Sensing, in press. 

The first was funded by a NASA grant issued to another PI and was submitted before the 
start of this grant period. However, it contains the full details of the statistical cloud 
analysis scheme and is included here for completeness. The second paper is an expansion 
of a presentation given at the IAMAP Symposium for Remote Sensing in the Polar 
Regions in Reading, England, in August 1989. It has been accepted for publication in a 
special issue of the International Journal of Remote Sensing devoted to the topic of remote 
sensing of the polar regions. 
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Analysis of Polar Clouds from Satellite Imagery Using Pattern Recognition 
and a Statistical Cloud Analysis Scheme 

Elizabeth E. Ebert 

Department o f Meteorology, University of Wisconsin, Madison . Wisconsin 
(Manuscript received 18 June 1988. in final form 26 November 1988) 

ABSTRACT 

The analysis of cloud cover in the polar regions from satellite data is more difficult than at other latitudes 
because the visible and thermal contrasts between the cloud cover and the underlying surface are frequently 
quite small. Pattern recognition has proven to be a useful tool in detecting and identifying several cloud types 
over snow and ice. Here a pattern recognition algorithm is combined with a hybrid histogram-spatial coherence 
(HHSC) scheme to derive cloud classification and fractional coverage, surface and cloud visible albedos and 
infrared brightness temperatures from multispectraJ AVHRR satellite imagery. The accuracy of the cloud fraction 
estimates were between 0.05 and 0.26. based on the mean absolute difference between the automated and 
manual nephanalyses of nearly 1000 training samples. The HHSC demonstrated greater accuracy at estimating 
cloud fraction than three different threshold methods. An important result is that the prior classification of a 
sample may significantly improve the accuracy of the analysis of cloud fraction, albedos and brightness tem- 
peratures over that of an unclassified sample. 

The algorithm is demonstrated for a set of AVHRR imagery from the summertime Arctic. The automated 
classification and analysis are in good agreement with manual interpretation of the satellite imagery and with 
surface observations. 


1. Introduction 

Understanding the nature and distribution of global 
cloud cover is a key element in improving our ability 
to understand and model the global climate. Indeed, 
the improv ed knowledge of global cloud fraction, height 
and thickness has been identified as a critical factor for 
improving weather and climate forecasts using general 
circulation models (WMO 1978: ECMWF 1981). To 
this end. the first project of the World Climate Research 
Programme is the International Satellite Cloud Cli- 
matology Project (ISCCP) (WMO 1982). which began 
collecting global satellite radiance data in 1983 and 
will continue until 1990 I Rosso w personal commu- 
nication). The goal of the ISCCP is to assemble a set 
of global radiance data over an interval of several years, 
from which 3-hourly cloud cover statistics may be de- 
rived and made available to the modeling community. 
One phase of the ISCCP was the development of a 
general procedure for estimating cloud cover from vis- 
ible and infrared radiances: this procedure is detailed 
by Rossow et ai. ( 1983). 

In the polar regions the climate is especially sensitive 
to the presence of clouds because of their strong radia- 
tive influence on the surface energy balance. The anal- 
ysis of cloudiness in polar high latitudes presents some 


( 'orresporuimg auinor address Elizabeth E. Ebert. Bureau of Me- 
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potential problems for the ISCCP algorithm, as well as 
many other cloud algorithms developed for use in lower 
latitudes. The absence of sunlight during winter renders 
the visible radiances useless during much of the year. 
During the summer the albedos of the clouds u.:d un- 
derlying snow and ice surface are similar, again re- 
stricting the usefulness of the visible information. The 
infrared radiances contain useful information during 
all seasons, but the distinction between cloud and sur- 
face is often unclear due to the frequent isothermal or 
inversion structure of the tropospheric temperature 
profile. Even the definition of cloudy and clear skies 
can be ambiguous. Cold season ice crystal precipitation 
(“diamond dust") from clear skies also have strong ra- 
diative effects in the lower troposphere (Curry 1987). 
and might justifiably be considered as cloudiness. 

Because the geostationary satellites do not ade- 
quately view the polar regions, polar orbiting satellites 
arc the only source of remotely sensed data for these 
regions. These data are less convenient to use than the 
geostationary' satellite data because of the moving frame 
of reference and reduced tem poral coverage, but this 
is offset by the increased spa' ui coverage and. on the 
TIROS-N and NOAA sateiln . the increased spectral 
resolution of the Advanced \ ary High Resolution Ra- 
diometer (AVHRR). Results from an ISCCP workshop 
on polar clouds (WMO 1987) suggest that the addi- 
tional spectral channels, beyond the visible and infrared 
channels common with the geostationary satellites, do 


c 1989 American Meteorological Society 


May 1989 


ELIZABETH E. EBERT 


7 S3 


indeed allow clouds to be distinguished from snow and 
ice. In particular, information from the near-infrared 
portion of the spectrum can be used to differentiate 
clouds of liquid phase from clouds of ice phase and 
the snow and ice surface ( Knottenberg and Raschke 
1982; Kidder and Wu 1984; Raschke et al. 1986; 
Molders 1987). and to detect low clouds over land or 
ocean at night (Eyre et al. 1984; d'Entremont and 
Thomason 1987). Solar radiation at 3.7 p m (AVHRR 
channel 3) is reflected well by cloud water droplets but 
not by the larger ice crystals found in cirrus clouds or 
by the surface (Arking and Childs 1985). 

The Textural appearance of the imagery can be used 
to help identify surface and cloud types. Cirrus, cu- 
mulus. and stratocumulus clouds have significantly dif- 
ferent macro- and microtextural characteristics which 
can be quantitatively measured (Welch et al. 1988; Kuo 
et al. 1988). Features such as cold, bumpy cloud tops, 
cloud shadows on the snow, illumination of cloud sides, 
and cracks and leads in the sea ice can aid the human 
observer in distinguishing cloud cover from sea ice and 
snow in satellite imagery (Kukla 1984; Welch et al. 
1989). Ebert ( 1987) used several spectral and textural 
features in an automated pattern recognition algorithm 
to classify 18 surface and cloud types with a skill of 
approximately 84T from Arctic summertime AVHRR 
imagery. Pattern recognition is yet a relatively new ap- 
proach to analyzing meteorological satellite imagery , 
and has been used successfully to measure cloud cover 
in other pans of the globe (e.g.. Hams and Barrett 1978; 
Bunting and Fournier 1980; Wu et al. 1985: Burhend 
et al. 1987; Seze and Desbois 1987: Garand 1988). 

It is likely that the a prion knowledge of a cloud 
type can provide information which would aid in the 
analysis of its radiative properties, and possibly also 
the radiative properties of the underlying surface. In 
this study the pattern recognition aigomhm of Ebert 
( 1987) is used to classify high latitude surface and cloud 
types from AVHRR imagery. The classification is then 
input to a hybrid histogram-spatial coherence analysis 
scheme which determines the cloud fraction, surface 
properties, and cloud albedo and bnghtness tempera- 
ture. Compansons against threshold methods with and 
without classification reveal the extent to which clas- 
sification can improve the cloud analysis. The com- 
bined algorithm is demonstrated for one day of sum- 
mertime satellite imagery from the Arctic. This data is 
pan of a Polar Cloud Pilot Data Set being analyzed 
simultaneously by several research groups in an effort 
to develop an improved scheme for detecting cloud in 
the polar regions for the 1SCCP (W MO 1987). 

The satellite data and auxiliary data for verifying the 
unalvses are described in section 2. Section 3 describes 
the combined classification and analvsis algorithm, as 
weil as its ability to correctlv diagnose the cloud type 
and traction lor a set of training samples. The classi- 
fication and analvsis of imagery from the Polar Cloud 
Pilot Data Set are presented and discussed in section 


4. The sensitivity of the results to underlying assump- 
tions in the algorithm, the propagation of error, and 
comparisons to threshold methods are made in section 

5, Section 6 presents a summary and discussion. 

2. Data 

a . Satellite data 

The algorithm was trained and tested using NOAA- 
7 AVHRR satellite imagery collected over the polar 
regions at the lower resolution Global Area Coverage 
(approximately 3 X 5 km held of view at nadir). The 
spectral data include AVHRR channels 1 (0.55-0.68 
^m), 2 (0.725- 1 . 1 /zm), 3 (3.55-3.93 and 4 ( 10.5- 
11.5 Mm), representing visible, near-infrared and in- 
frared radiances. [The data from AVHRR channel 5 
(11.5 ^um-12.5 were not included because the 
AVHRR instruments on many polar orbitting satellites 
(TIROS-N. NOAA-6. and NOAA-8) do not measure 
radiances in this spectral band.] These data are read 
and processed from tape using the Man-computer In- 
teractive Data Access System (McIDAS) at the Uni- 
versity of Wisconsin. Calibration of the pixel bright- 
nesses to albedos, p (channels l. 2. and 3). and bright- 
ness temperature T B (channel 4) is made according to 
Kidw'ell (1984). The albedo in channel 3. which con- 
tains both reflected solar and emitted thermal radiation, 
is estimated by subtracting the thermal radiance that 
would be emitted between 3.55-3.93 ^m from a black- 
body radiating at the brightness temperature measured 
in channel 4. Since the composition of the scene is no* 
known a priori, this is only a rough approximation 
the true channel 3 albedo. It will tend to slightly un- 
derestimate the 3.7 pm albedo of overcase low clouds, 
but overestimate the 3.7 albedo of broken cirrm 
clouds. In general, however, these errors are not larz 
enougn to cause confusion between water and 
clouds. 

Each image is divided into an array of cells of di- 
mension 2.5° latitude by 2.5° longitude, in accordance 
with the 1SCCP Polar doud Pilot Study ( WMO 1987). 
Near t>0°N, this comprises about 2400 pixels, while at 
the pole a 2.5° X 2.5° cell has only 1 1 2 pixels. In order 
to make the cell areas more nearly equal and provide 
enough pixels (i.e.. data points) for the classification 
and histogram analysis, the ceils nearest the pole are 
combined for those steps. The cloud fraction, however, 
is still analyzed for 2.5° n 2.5° cells everywhere.) Be- 
tween 80°N and 85 °N the cell width increases to 5 3 
longitude, between 85 C N and 87.5°N the cell width 
again doubles to 10° longitude, and between 87.5°N 
and the pole the cell width again doubles to 20° lon- 
gitude. The 45 pixels at the extreme left and right of 
each satellite image are discarded because of significant 
distortion due to the earth's curvature. 

The temporal coverage of the satellite imagery 
spanned the dates from I — July |984. Complete cov- 
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Table 1 . Auxiliary data for verifying the Arctic cloud analyses. 

Data type 

Location 

Source 

Surface synoptic observations* 
Surface synoptic observations 
Surface synoptic observations 
Hemispheric analyses of sea level 
pressure and 500 mb 
geopotential height 
Weekly sea ice concentration 

North American Arctic 
European Arctic 

Fram Strait (approx. 80°N, 0°E) 
Northern Hemisphere 

Northern Hemisphere 

National Atmospheric and Oceanic Administration 
Berliner Wetterkarte (publ. Free University of Berlin) 
Lindsay (1985) 

European Centre for Medium Range Weather 
F- recasting 

Navy-NOAA Joint Ice Center, Suitland, Maryland 


* Surface synoptic observations include sea level pressure, surface temperature and dewpoint, wind velocity, surface weather, cloud 
fraction, and frequently cloud type. 


erage of the Arctic from 60° to 90° latitude and 90°W 
to 90°E longitude, plus additional incomplete coverage 
between 90°E and 90°W, were included. A subset of 
this data representing a variety of surface and cloud 
types was selected for training and testing the algorithm. 
The training sample dimensions were 32 X 32 pixels, 
which were chosen to simplify the testing of certain 
textural features while providing an area compatible 
with that used for the cloud analysis. 

b. Auxiliary data 

While the classification and analysis algorithm em- 
ploys satellite data alone, supplementary data are 
needed to verify the algorithm's analyses. There is no 
extensive network of surface observing stations in the 
Arctic, so additional comparisons must be made against 
the scattered data which are available, including surface 
observations from isolated stations. Navy sea ice esti- 
mates. and special observations made during the 1984 
Marginal Ice Zone Experiment {MIZEX 84). Table 1 
summarizes the auxiliary data used to verify the cloud 
analyses. 

3. Cloud analysis algorithm 

The algorithm comprises two steps for each 2.5° 
X 2.5° box. First, a classification algorithm identifies 
the surface and cloud type within the box. This infor- 
mation is then used in a hybrid histogram-spatial co- 
herence analysis procedure which estimates cloud 
fraction, and clear and cloudy albedos and brightness 
temperatures. 

a. Classification of the scene 

The classification algorithm is based on the pattern 
recognition algorithm described by Ebert ( 1987) and 
shown schematically in Fig. 1. A summary of this clas- 
sification algorithm is presented here. 

Spectral and textural features characterizing each cell 
are extracted from the radiance data. Spectral features 
describe the characteristics of the sample as a whole, 
such as its mean albedo or its brightness temperature 
range. The textural features describe aspects of the 


variability, or “bumpiness,” of a scene, and are com- 
puted from the albedo or brightness temperature dif- 
ferences between each individual pixel and its neigh- 
bors. 

The classification algorithm employs the following 
numerical features, in order of importance: 1) Mean 
channel 1 albedo, 2) Mean channel 4 brightness tem- 
perature, 3) Channel 1 albedo - channel 2 albedo. 4) 
Percentage of pixels with < 0.08, 5) Percentage of 
pixels with T B > 273, 6) Ratio of channel 3 and channel 
1 albedos, 7) Percentage of pixels with p, < 0.15, 8) 
Mean channel 4 angular second moment. 9) Maximu-' 
channel 3 entropy, 10) Maximum channel 1 entrop>. 
At least one spectral and one textural feature were cho- 
sen in each channel. 1 The first seven features desenbe 
the warm, dark characteristics of land and water sur- 
faces, and the cooler, brighter nature of the clouds and 
the snow and ice surface. The fourth and sixth features, 
which make use of the channel 3 albedo, aid in distin- 
guishing the highly reflective water-phase clouds from 
the poorly reflecting ice crystal clouds and the snow 
and ice surface. The last three features are textural fea- 
tures, and generally discriminate between smooth and 
variable fields. 

The classification of a sample is made using the 
maximum likelihood decision rule, a statistical tech- 
nique which is widely employed in classification prob- 
lems. The maximum likelihood decision rule assumes 

the existence of distinct classes. C, , C 2 Cj, which 

may be represented by clusters in the 10-dimensional 
feature space. In essence, a sample or observation which 
is characterized by a i 0-dimensional feature vector, v 

= [tv v : where each feature. i\, describes 

some aspect of the sample (mean channel l albedo, 
etc.), will be assigned to that class, C k , to which its 
feature vector is closest in the feature space. 

Eighteen surface and cloud categories were differ- 
entiated in the training data: l ) Land without snow. 
2) Land with snow. 3) Open water. 4) Unbroken sea 


1 This feature list differs slightly from that of the Ebert < 1987) al- 
gorithm due to improved calibration procedures used in the more 
recent version. 
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Multispectral AVHRR imagery 



Fig. 1. Flowchart of the automated classification algorithm. 


ice. 5) Broken sea ice. 6) Ice sheet. 7) Mixed surface 
types. 8) Cs/As over land/water, 9) Ci/Ci/Ac over land/ 
water. 10) St over land/water. 1 1) Sc over land/water, 
12) Cu over land/water. 13) Cs/As over snow/ice. 14) 
Ci/Cc/Ac over snow/ice. 1 5) St over snow/ice. 16) Sc 
over snow/ice. 17) Cu over snow/ice. 18) High cloud 
over low cloud. This selection of categories is somewhat 
subjective, and different investigators might opt for dif- 


ferent categories. In fact, many classification schemes 
use unsupervised clustering techniques to allow the data 
to fall into “natural” categories. The danger there, 
however, is that the classes may not necessarily rep- 
resent logical physical groupings, and the clustering is 
highly dependent on the distribution of samples. 

Here the classes are specified as certain surface and 
cloud types which have geophysical significance. The 
clouds are separated according to whether they are 
stratiform or broken, ice or water phase, and are over 
a land/ocean or snow/ice surface. The distinction be- 
tween surface types underlying the cloud classes is made 
because the differing visible and thermal properties of 
the warm, dark land and water surfaces, and cold, 
bright snow and ice can influence both the albedo of 
partially transparent clouds due to multiple reflections, 
and the cloud brightness temperature due to trans- 
mission of thermal radiation from the surface. Cirri- 
form and alto clouds are defined here to be composed 
primarily of ice crystals, and stratus, stratocumulus and 
cumulus are defined to be composed primarily of liquid 
water droplets. Class 17, cumulus over snow/ice, refers 
not only to convective cumuliform clouds, but also to 
small (but not subgrid scale) isolated clouds which may 
be shallow in depth. Class 18 refers to clouds of ice 
phase overlying clouds of liquid phase. 

An iterative training procedure was used to teach 
the algorithm the characteristics of the 1 8 classes, based 
on the features derived from the training samples. The 
algorithm's skill, 5, at classifying the (training samples 
was 85%, as defined by ! 

•S' = ~~~ X 100 (%), (1) 

T — t 

where R is the number of correct classifications. T is 
the total number of classifications, and E is the number 
of correct classifications expected due purely to random 
chance. Aside from layered cloudiness, the most com- 
mon classes present in the training imagery were stratus 
and stratocumulus over land/water and snow/ice, typ- 
ical for the Arctic during summer (Curry and Herman 
1985). The algorithm had the most difficulty identifying 
high cloud over low cloud, misclassifying it as either 
low cloud or high cloud alone. Cirrus over land/water 
was sometimes undetected, as was stratus cloud over 
sea ice on some occasions. The algorithm also failed 
to find very small amounts of cumulus cloudiness over 
land. Overall, about 2% of the cloudy regions were 
misclassified as clear, and less than 1% of the clear 
regions were misclassified as being cloudy. The esti- 
mated skill of the algorithm in classifying independent 
imagery is comparable to that of the training set when 
verified against the analyses of independent experts 
(Ebert 1987). 

b. Analysis of the scene 

Information from the classification phase is com- 
bined with a hybrid histogram-spatial coherence 
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FlG. 2. Schematic diagram of the hybrid histogram-spatial 
coherence method. 


method adapted from Liu et al. (1988) to determine 
the fractional cloudiness of the scene. This scheme is 
ideally suited for analyzing ensembles of pixels. The 
cloud fraction. N, is defined here to represent the frac- 
tion of an area covered by cloud when viewed from 
directly above. For clear classes with two types of sur- 
faces, such as snow on land, broken sea ice, and mixed 
surface types, the fractional amount of each surface 
type is computed. The channel 1 albedo and channel 
4 brightness temperature of the clear and cloudy por- 
tions of the scene (denoted T v , a c , and T t ) are de- 

termined as part of the analysis. These represent sat- 
ellite observed values (i.e., uncorrected for atmospheric 
effects), which may differ from the true surface and 
cloud albedos and temperatures. 


1) Determination of cloud fraction 

The hybrid histogram-spatial coherence method 
(HHSC) is shown schematically in Fig. 2. The object 
of the scheme is to identify clear, partly cloudy, and 
completely cloudy pixels within a scene and to deter- 
mine the mean albedos or brightness temperatures as- 
sociated with the completely clear and completely 
cloudy pixels. The cloud fraction is then easily obtained 
using the spectral channel which shows the greatest 
contrast between cloud and surface. 

The optimal channel for determining cloud (ice) 
fraction for each class is given in Table 2. Each of the 
AVHRR channels used was optimal for at least one 
class. Channel 2 was employed for mixed surface types 
because land, water, and snow/ice each have a distinc- 
tive signature in this spectral region. For “high cloud 
over low cloud,” channel 3 is used to estimate the low 
cloud fraction, and channel 4 is used to determine the 
high cloud fraction. 

First we use a histogram technique in the class-op- 
timal channel to find the albedos or brightness tem- 
peratures representing completely clear and completely 
cloudy pixels in the scene. The histogram technique is 
based on the assumption that, due to random instru- 
ment noise and natural variance, the values of p or T B 
in an ensemble of clear or cloud filled pixels are dis- 
tributed normally, with the peak of the distribution 
representing the mean value. (We note that, except for 
channel 3, the AVHRR data is nearly noise-free, so 
that the observed variation is due almost entirely to 
natural variance. In this case the distribution of the 
pixel values is not necessarily Gaussian; nevertheless, 
we find that, in practice, this technique can quite ad- 
equately determine a representative value of p or T B 
corresponding to the ensemble of pixels.) For a sample 
containing both clear and cloudy regions, we expect to 
find two peaks in the histogram, corresponding to the 
clear and cloudy albedos or brightness temperatures. 
The use of the class-optimal channel in this step max- 
imizes the separation of these peaks. Partly cloudy pix- 
els obscure the shape of the distribution curve between 
the clear and cloudy peaks. Therefore the geometry of 
the outer portions of the distnbution is used to infer 


Table 2. Class optimal AVHRR channel for Arctic summer cloud fraction analysis. 


Class Optimal channel Class Optimal channel 


1. Land l 

2. Snow on land 1 

3. Open water 1 

4. Unbroken sea ice 1 

5. Broken sea ice 1 

6. Ice sheet 1 

7. Mixed surface types 2 

8. Cs/As over land/water 4 

9. Ci/Cc/Ac over land/water 4 


10. St over land/water 1 

1 1 . Sc over land/water l 

12. Cu over land/water I 

13. Cs/As over snow/ice 4 

14. Ci/Cc/Ac over snow/ice 4 

1 5. St over snow/ice 3 

16. Sc over snow/ice 3 

1 7. Cu over snow/ice 3 

1 8. High cloud over low cloud 3. 4 
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the means of the clear and cloudy peaks. In the visible 
channel this would correspond to the dark side of the 
surface albedo distribution and the bright side of the 
cloud albedo distribution. 

Each peak in the histogram is analyzed using the 3- 
point method (Brower et al. 1976). For a normal prob- 
ability density function for albedo or brightness tem- 
perature, x y the frequency of occurrence of x is given 
by 

fix) = fo exp{-(.v - fid ) 2 /2a 2 }, (2) 

where the peak frequency, mean and standard devia- 
tion {fo, mo and a, respectively) are unknown. By taking 
three points (x, ,f) and (x k ,f k ) from the outside 

of the distribution curve (see Fig. 2), three equations 
in three unknowns. 

f = Jo exp{-(.v, - mo) 2 /2 a 2 } 

fj=fo exp{— (x) - Mo) : / 2 <r} 

fk =fo exp {-(•** - no ) 2 /2a 2 }, (3) 

can be solved using Gaussian elimination to obtain the 
peak value 

= XJ_ log UfJk) - xj log Wfk) + X% log Wf) 

^ 2[X, log (fj/fk) - Xj log Wfk) + x k log if If,)) 

and the standard deviation 

= ~[U- ~ x) 2 - Imix, - .v,)] 

* 2 log (M) 


For each peak, values of and <r are calculated from 
many different 3-point combinations, then plotted 
against each other on a two-dimensional histogram. 
The combination of (jio, <*) with the greatest frequency 
of occurrence in the 2-D histogram is taken to be rep- 
resentative of the mean clear (cloudy) value, il s (ix c ), 
and its standard deviation, a s ( a c ), for the cell. 

It is necessary to specify reasonable characteristic 
values and bounds on pixel albedos and brightness 
temperatures in order to avoid misidentification of clear 
or cloudy peaks. These bounds are initially estimated 
from the characteristic mean, maximum and mini- 
mum values for the classes in the training set, and con- 
tinuously updated during the analysis procedure. Initial 
class/channel characteristic values for the cloud analysis 
algorithm are given in Table 3, along with bounds 
(thresholds) on pixel values which can be considered 
completely clear or completely cloudy. 

If only one peak is clearly defined in the histogram, 
as would be the case for a cloud fraction near 0.0 or 
1 .0, then its definition as either “surface” or “cloud” 
is dependent on whether it more nearly matches the 
characteristic clear or cloudy value for that channel 
and class. In that case, the most likely type (surface or 
cloud) is assigned the peak value, while the other type 
is assigned its characteristic value for the class. If no 
peaks are definable, as for clouds of varying tempera- 
ture or optical thickness covering the scene, the class 
characteristic values for both the surface and cloud 
values are used. 

Next, the spatial coherence method (Coaklev and 
Bretherton 1982) is used to isolate completely clear 


Table 3. Initial class characteristic values for surface and cloud for the Arctic summer cloud analysis algorithm. The values in parentheses 
are the bounds on clear and cloudy albedos and brightness temperatures which separate the two. For classes 2, 5. and 7, the “cloud’' values 
refer to a second surface type, generally snow/ice. 


1. Land 

2. Snow on land 

3. Open water 

4. Unbroken sea ice 

5. Broken sea ice 

6. Ice sheet 

7. Mixed surface types* 

8. Cs/As over land/water 

9. Ci/Cc/Ac over land/water 

10. St over land/ water 

1 1 . Sc over land/water 

12. Cu over land/water 

13. Cs/As over snow/ice 

14. Ci/Cc/Ac over snow/ice 

15. St over snow/ ice 

16. Sc over snow/ice 

17. Cu over snow nee 

18. High cloud over low cloud (low) 

(high) 


Ch. 1 

Clear Cloud 

(albedo X 100) 

Ch. 3 

Clear Cloud 

(albedo x 100) 

Ch. 4 

Gear Cloud 

(K) 

11 (17) 




285 (293) 


II (17) 

64 (54) 



285 (276) 

273 (276) 

8 (15) 




279 (287) 


49 (64) 




271 (274) 


8 (15) 

49 (32) 



279 (271) 

271 (275) 

64 (74) 




259 (270) 


10 (17) 

56 (32) 

11 (21) 

27 (10) 

282 (271) 

265 (275) 

10 (17) 

56 (31) 



282 (271) 

248 (275) 

10 (17) 

45 (21) 



282 (271) 

248 (275) 

10(17) 

45 (26) 



282 (271) 

269 (282) 

10(17) 

55 (20) 



282 (271) 

266 (282) 

10(17) 

53 (20) 



282 (271) 

265 (282) 

56 (74) 

69 (35) 



265 (255) 

251 (270) 

56 (74) 

69 (22) 



265 (255) 

251 (273) 

56 (74) 

61 (39) 

2 (5) 

17 (5) 

265 (255) 

262 (274) 

56 (74) 

71 (31) 

2 (5) 

28 (10) 

265 (255) 

255 (279) 

56 (74) 

71 (31) 

2 (5) 

23 (10) 

265 (255) 

255 (273) 

56 (74)** 

55 (21) 

2 (5) 

22 (5) 

265 (255)** 

263 (282) 


60 (21) 




252 (275) 


* For class 7. the class characteristic values listed under Ch. 3 are actually for Ch. 2. 

** Equatorward of 75°N. the surface characteristic values are initialized to those of classes 8-12. 
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and completely cloudy pixels. This corresponds to Liu 
et al.’s (1988) “improved two-threshold method.” Local 
means and standard deviations of p or Ty in 4-pixel 
squares within a cell are plotted against each other to 
form an arch-like structure. The feet of the arch rep- 
resent completely clear and completely cloudy regions, 
while the central portion of the arch contains partly 
cloudy pixels. An arch “foot” is defined here as the 
region containing local mean values within two stan- 
dard deviations (2 a) of no as determined from the his- 
togram procedure, and whose local standard deviations 
are also less than 2a. In cases where the feet overlap, 
a threshold is set which separates the two regions in 
proportion to their relative values of a. In addition, 
the thresholds separating completely clear and com- 
pletely cloudy pixels (i.e., p±2<r) were constrained not 
to go beyond their respective bounds given in Table 3. 

Pixels identified as completely clear or completely 
cloudy are tagged and given weights ot w = 0 or w = 1 , 
respectively, where the weights represent pixel cloud 
fraction. Those pixels outside the arch feet are given 
weights 


YV — , v “ ' 

Me - Mi 

subject to the constraint that 0 < w < l . The weights 
are then summed over the entire scene and divided by 
the total number of 4-pixel squares to obtain the cloud 
fraction. N. 

Although the selection of 2 a to define the arch teet 
may seem rather lax, it is necessary to guarantee enough 
pixels for the subsequent histogram analysis of the re- 
maining surface and cloud properties, described in the 
next subsection. Moreover, we have found that the de- 
rived cloud fraction is not very sensitive to the choice 
of 2a when compared to the results obtained using 1 a 
or 0 (i.e., no arch feet). 

2) Determination of remaining surface and 

CLOUD PROPERTIES 

Once the completely clear and completely cloudy 
pixels have been tagged using the HHSC in the class- 
optimal channel, these pixels are used to generate 
"clean” histograms in the remaining spectral channels. 
The three-point histogram analysis is performed on 
the "clean” histograms to obtain a v , a t , T s and T c . 
There is an advantage to combining the histogram 
method with the spatial coherence method, rather than 
using the histogram method alone or the spatial co- 
herence method alone in each spectral channel. The 
natural variability ot the GAC pixels over the 2.5 
X 2.5° region causes the teet ot the spatial coherence 
arch to be poorly defined. The histogram method is 
able to find mean values for the clear and cloudy peaks 
in most cases. This is used to define the regions ot the 
feet for the spatial coherence method, which sorts the 
clear and cloudy pixels. This sorting better defines the 


peaks in the remaining channels, yielding a more ac- 
curate estimate of the mean surface and cloud prop- 
erties in those channels. 

Often one peak dominates the histogram. This re- 
sults when the clear and cloudy values are nearly equal, 
as for clouds overlying snow and ice, or when only a 
small amount of surface or cloud is present. To resolve 
this problem, the pixels in the "clean” histogram are 
split according to the cloud fraction determined in the 
optimal channel, and each side of the histogram is an- 
alyzed independently. In this way, a minor peak rep- 
resenting a small number of cloudy pixels will not be 
ignored in favor of a larger surface peak. Again, if a 
clear or cloudy peak cannot be found, its class char- 
acteristic value is used. 

The characteristic values are updated after each cell 
is analyzed by keeping a running mean of clear and 
cloudy values for each channel and class. This ensures 
that the characteristic values are representative of the 
surface and cloud values already encountered in the 
scene. Only at the beginning of each satellite image are 
these characteristic values reinitialized to prespecified 
values determined from the training set. 

c. Testing the cloud analysis algorithm 

An objective test of the HHSC was conducted by 
creating and analyzing synthetic cloud samples. These 
samples were generated from the 4-channel class char- 
acteristic means and standard deviations using a Gaus- 
sian random number generator, and arranged in the 
six test patterns shown in Fig. 3: (a) checkerboard, (b) 
overcast, (c) cloud edge, (d) complete gradient, (e) par- 
tial gradient, and (0 sine wave. The test patterns ap- 
proximate cloudy situations which occur in the im- 
agery. Ten simulations were run for each cloud type 
and test pattern. 

Table 4 shows the means of the simulations for three 
cloud classes, each with a different class-optimal chan- 
nel. The analyzed cloud fraction is generally correct to 
within 0. 10, and always within 0.20. The analyzed sur- 
face and cloud albedos and brightness temperatures 
are usually accurate to within about 0.05 and 4 K. 
respectively. The most difficult scenes to analyze ac- 
curately were the overcast and complete gradient pat- 
terns. The underestimation ot the overcast cloud trac- 
tion occurred because the pixels at the dark or warm 
end of the histogram were analyzed as being partly 
cloudv instead of completely cloudy. The scene com- 
posed entirely of a horizontal gradient, (d), was prob- 
lematic because the histogram contains no obvious 
surface or cloud peaks. The clear and cloudy albedos 
and brightness temperatures analyzed by the algorithm 
tend to be closer to the center ot the histogram than 
the true values are. The overall pertormance ot the 
HHSC on the artificial cloud fields is summarized in 
Table 5 for all of the classes. Except for the second test, 
the algorithm slightly overestimated cloud fraction, 
typically by about 0.05-0.10. 
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Fig. 3. Test patterns and associated histograms for six 
synthetic cloud samples. 


The HHSC was also tested against a manual neph- 
analysis of the training samples. Figure 4 shows the 
manual and automated analyses of cloud fraction for 
the three training images. (These correspond to Fig. 2 


of Ebert 1987.) There is generally good qualitative 
agreement between the two analyses. The clear and 
cloudy regions are well diagnosed, especially in the sec- 
ond and third images. Table 6a quantitatively compares 
the manual and automated analyses for all of the cloud 
samples, including those which were misclassified. As- 
suming the manual analyses to be accurate, the algo- 
rithm underestimated fractional cloud cover for five 
of the classes, and overestimated it for six classes. The 
largest errors occurred in the analysis of high cloud 
fraction as measured by the root-mean-square (RMS) 
error, and the cumulus cloud fractions were analyzed 
most accurately. The mean absolute error averaged 
about 0.13 for low clouds and 0.23 for high clouds. 

By comparing the manual nephanalysis to the au- 
tomated nephanalysis for only those cloud samples 
which were correctly classified, one can determine the 
accuracy of the cloud algorithm under the best of cir- 
cumstances. The error statistics for the correctly clas- 
sified samples are presented in Table 6b. The RMS 
errors do not differ greatly from those obtained using 
all the samples. Less than 3% of the cloudy samples 
were misclassified as clear, while the remaining mis- 
classifications were between various cloud types. The 
HHSC does not appear to be extremely sensitive to 
these misclassifications. (A more thorough discussion 
of algorithm sensitivity to classification errors appears 
in section 5.) 

The tendency of the hybrid histogram-spatial co- 
herence method to predict low cloud fraction more 
accurately than high cloud fraction is in contrast to 
the performance of the U.S. Air Force three-dimen- 
sional nephanalysis (3DNEPH), which relies on sat- 
ellite, surface and aircraft observations (Fye 1978). In 
a comparison of the 3DNEPH cloud fraction estimates 
to manual estimates using high resolution DMSP im- 
agery over the Beaufort Sea in June 1980, Curry and 
Herman (1985) found a mean absolute difference of 
0.43 for low clouds, 0.12 for middle clouds, and 0.06 
for high clouds. This indicates the difficulties associated 
with using an infrared threshold technique to detect 
low clouds at polar latitudes. 


4. Analysis of Polar Cloud Pilot Data 

In this section, the surface and cloud cover are an- 
alyzed and discussed for l July 1984, the first day of 
the Arctic summer imagery in the Polar Cloud Pilot 
Data Set. Fourteen satellite images cover the region 
between 60° and 90° north latitude. At least one sat- 
ellite overpass was used for all regions of the Arctic 
except the Alaskan and far eastern portions, and a small 
region immediately surrounding the pole, where the 
pixels were too close to the edge of the image and were 
discarded due to distortion (see section 2). Each image 
was analyzed individually, then the analyses combined 
and averaged to form a hemispheric mosaic. 
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Table 4. Results of the cloud analysis simulations for three cloud classes, each with a different class optimal channel: a, = ch. 1 surface 
albedo; p 3j - ch. 3 surface albedo; T s = surface brightness temperature; a c = ch. 1 cloudy albedo; p* = ch. 3 cloudy albedo; T c = cloudy 
brightness temperature. 


(a) Class 1 1: 

: Sc over land/water (a, = 0.10; a c = 

0.55; T s 

= 282; T c = 

266; Class optimal channel 

= i) 


Test 

Actual N 

Algorithm N 

ct s 



T, 

T c 

(a) Checkerboard 

0.55 

0.45 


0.09 

0.60 


282 

271 

(b) Overcast 

1.00 

0.88 


0.10 

0.61 


282 

264 

(c) Cloud edge 

0.15 

0.14 


0.09 

0.55 


282 

266 

(d) Complete gradient 

0.52 

0.62 


0.14 

0.43 


279 

268 

(e) Partial gradient 

0.36 

0.36 


0.09 

0.56 


282 

267 

(f) Sine wave 

0.58 

0.59 


0.09 

0.51 


281 

265 

(b) Class 15: St over snow/ice (a, * 0.56; 

a c = 0.61; p3, = 0.015; p^ : 

= 0.165; T, = 

: 265; T c = 262; Gass optimal channel = 3) 


Test 

Actual N 

Algorithm N 

a. 


Pis 

Pic 

T, 

T c 

(a) Checkerboard 

0.50 

0.52 

0.56 

0.61 

0.014 

0.187 

265 

262 

(b) Overcast 

1.00 

0.86 

0.58 

0.63 

0.015 

0.182 

264 

261 

(c) Goud edge 

0.15 

0.21 

0.54 

0.61 

0.015 

0.166 

265 

262 

(d) Complete gradient 

0.52 

0.68 

0.55 

0.59 

0.032 

0.101 

265 

262 

(e) Partial gradient 

0.36 

0.47 

0.55 

0.61 

0.013 

0.148 

264 

261 

(0 Sine wave 

0.58 

0.63 

0.55 

0.60 

0.017 

0.138 

265 

262 

(c) Class 9: Ci/Cc/Ac over land/water (a, = 0. 10; a c 

= 0.45; T s = 282; T c 

= 259; Gass optimal channel = 4) 


Test 

Actual N 

Algorithm N 




T, 

T c 

(a) Checkerboard 

0.50 

0.51 


0.10 

0.47 


285 

256 

(b) Overcast 

1.00 

0.95 


0.10 

0.48 


281 

255 

(c) Goud edge 

0.15 

0.23 


0.10 

0.45 


285 

266 

(d) Complete gradient 

0.52 

0.66 


0.14 

0.37 


278 

267 

(e) Partial gradient 

0.36 

0.43 


0.10 

0.47 


285 

260 

(f) Sine wave 

0.58 

0.61 


0.11 

0.46 


283 

262 


Mosaics of visible and thermal infrared imagery for 
the 1 July 1984 Arctic are shown in Fig. 5. The imagery 
shows a large cyclone in the western Canadian Arctic 
and low cloud cover in the Canadian Archipelago and 
northwestern Greenland. The Greenland Sea contains 
broken sea ice, and further south a large region of stra- 
tus and stratocumulus blankets the ocean. Layered and 
cumulus clouds are present over Scandinavia, while a 
large region in Siberia is clear. The Arctic Basin con- 
tains multiyear sea ice with large areas of cloud cover. 

The image classification, presented in Fig. 6, and the 
analyzed cloud fraction, presented in Fig. 7, together 
describe the cloud cover analyzed by the algorithm. 
Overall, for the 1 July 1984 Arctic, the estimated frac- 


tional area covered by low clouds alone was 0.26, the 
estimated fractional area covered by high cloudiness 
alone was 0.05, and the fractional area covered by lay- 
ered cloudiness was estimated at 0.24. The most com- 
mon cloud types analyzed were stratus, both over water 
and over ice, and layered cloudiness. This matches with 
the findings of many previous studies ( Vowinckel 1 962: 
Gavrilova 1963; Huschke 1969). 

Many of the cloud features are related to the synoptic 
conditions, shown as analyses of sea level pressure and 
500 millibar geopotential height for the 1 July Arctic 
in Fig. 8. The comma cloud in the Northwest Terri- 
tories is part of an occluded cyclone. The line of cirrus 
clouds running to the northeast of Novaya Zemlya is 


Table 5. Summary of the results of the cloud analysis simulations. 


Test 

Number of 
samples 

Actual 

fraction 

Mean algorithm 
fraction 

RMS 

error 

Mean absolute 
difference 

(a) Checkerboard 

10 

0.50 

0.52 

0.07 

0.05 

(b) Overcast 

10 

1.00 

0.88 

0.13 

0.12 

(c) Cloud edge 

10 

0.15 

0.21 

0.09 

0.06 

(d) Complete gradient 

10 

0.52 

0.64 

0.13 

0.12 

(e) Partial gradient 

10 

0.36 

0.44 

0.10 

0.08 

(0 Sine wave 

10 

0.58 

0.61 

0.06 

0.05 








Fig. 4. Manual and automated cloud fraction analyses of the three Arctic summer training images. Cloud fraction is expressed in units of tenths. 
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T A m p 6 Statistics for the manual and automated nephanalyses for the Arctic summer cloudy training samples: 
TABLE 6 ‘ ^T^lmSL^lcorrecUy class, fled samples only. Refer to text for category numbers. 


(a) All training samples. 


Number of 


Mean manual 


Gass 

samples 

fraction 

g 

17 

0.85 

9 

34 

0.35 

10 

129 

0.79 

\ i 

60 

0.55 

12 

45 

0.21 

13 

12 

0.78 

14 

51 

0.49 

15 

98 

0.74 

16 

48 

0.76 

17 

14 

0.18 

18 

178 

0.79 


Mean algorithm 
fraction 


0.68 

0.47 

0.75 

0.60 

0.22 

0.83 

0.40 

0.63 

0.66 

0.21 

0.80 


RMS 

error 


0.34 

0.26 

0.15 

0.15 

0.11 

0.24 

0.35 

0.28 

0.18 

0.06 

0.23 


Mean absolute 
difference 


0.24 

0.22 

0.11 

0.11 

0.09 

0.18 

0.26 

0.19 

0.14 

0.05 

0.17 


(b) Correctly classified training samples. 

Number of 

Class samples 


8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 


13 

30 

114 

44 

34 

12 

49 

87 

44 

13 

130 


Mean manual 
fraction 

0.89 
0.35 
0.82 
0.55 
0.22 
0.78 
0.49 
0.77 
0.79 
0.18 
0.80 


Mean algorithm 
fraction 

0.66 
0.48 
0.88 
0.61 
0.25 
0.83 
0.40 
0.70 
0.67 
0.21 
0.88 


RMS 

error 

0.37 

0.27 

0.15 

0.14 

0.09 

0.24 

0.36 

0.24 

0.18 

0.06 

0.21 


Mean absolute 
difference 

0.26 
0.22 
0.11 
0.11 
0.07 
0.18 
0.27 
0.16 
0.14 
0.04 
0.15 


associated with the flow of the jet stream around an 
omega high, which produces subsidence and clearing 
near 90°E. In the analysis, the Norwegian Sea and 
North Atlantic are covered by dense stratus cloud cover 
(N > 0.80), with a small clear region over Iceland as- 
sociated with a surface high pressure center. Greenland 
is analyzed as mostly cloudy, and stratus and layered 
clouds are analyzed over much of the Canadian Ar- 
chipelago. Northern Siberia is analyzed as mostly clear, 
with scattered cumulus and cirrus clouds. The Arctic 
Ocean has primarily cloud-free pack ice in the Euro- 
pean sector, with stratus and layered cloudiness ana- 
lyzed in the Soviet and western regions. Much ot this 
cloudiness mav be associated with the southeasterly 
flow of warm moist air near 90° W (Jayaweera 1982; 

Curry and Herman 1985). 

The automated analysis can be compared to the sur- 
face and weather conditions observed on 1 July 1984. 
plotted in Fig. 9 for several stations. These observations 
consist of station reports from the Canadian and Eu- 
ropean Arctic, special observations Irom the M1ZEX 
study area, and the Navv-NOAA Joint Ice Center sea 
ice analysis from 3 July 1984. Before using these ob- 
servations to validate the satellite automated analyses, 
some limitations related to differences in temporal and 
spatial sampling must be addressed. The surface 


weather reports plotted in Fig. 9 were made at 3 UTC. 

6 UTC. 9 UTC and 15 UTC (depending on the data 
source) and may not necessarily represent conditions 
for the entire day, although the diurnal variation in 
surface and cloud conditions is small in the Arctic at 
this time of vear. Even comparisons between nearly 
concurrent measurements are impaired by the fact that 
the surface report is a point measurement, influenced 
by local topography and a limited sky view, while the 
satellite analysis covers a 2.5° X 2.5° region. In addi- 
tion. most stations in the polar regions are located on 
coastlines, transition zones between maritime, conti- 
nental. pack ice. and ice sheet regimes, and may thus 
experience unique weather conditions. Considering 
these limitations there appears to be good agreement 
between the analyzed and observed surface and cloud 
types in the 1 Julv Arctic. In particular, the cirrus 
cloudiness near 90°E is well analyzed, as are the layered 
clouds over Scandinavia and the partly cloudy condi- 
tions over much of the Canadian Archipelago. The an- 
alyzed cloud cover over Greenland and the pack ice 
cannot" be verified due to the lack ot surface observa- 
tions. . . 

The correct diagnosis of surface type is important 
in correctly estimating the surface radiation budget. 
This is particularly true during summer, as the amount 
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Fig. 6. Automated classification of the Arctic 1 July 1984 imagery. 


for this purpose. Figure 10 shows the analyzed sea ice 
cover for the 1 July Arctic. The automated analysis 
shows a small amount of sea ice in the Kara Sea, while 
west of Novaya Zemlya the ocean is analyzed as ice 
free. The sea ice edge curves over the northern side of 
Svalbard and becomes quite broken in the Greenland 
Sea. Broken sea ice is also found in Baffin Bay. Un- 
broken sea ice is analyzed throughout the polar basin. 

The algorithm’s analyzed sea ice fraction was verified 
against the Navy-NOAA sea ice concentration esti- 
mated for each 2.5° X 2.5° box from the 3 July ice 
chart. The expected change in sea ice extent during the 
two days separating the 1 July Polar Cloud Pilot Data 
and the Navy-NOAA observations is small, not ex- 
ceeding 15 km. When compared to the Navy-NOAA 
analysis, the algorithm accurately diagnosed the sea ice 


Fig. 5. Mosaic of visible (a) and infrared fb) imagery over the 
Arctic on l July 1984. fit should be noted that /or display purposes 
only the mosaics of the visible and thermal imagery were constructed 
using a maximum brightness criterion, and thus appear colder and 
cloudier than the individual images, due primarily to cloud motion 
between consecutive satellite overpasses.) 


of absorbed solar radiation is much larger for water 
and land, which have low albedos, than for the highly 
reflective snow and ice surface. The sea ice concentra- 
tion and position of the ice margin can be estimated 
by combining the results from the classification and 
analysis steps. Class 4 (unbroken sea ice), class 6 (ice 
sheet), and classes 13-17 (clouds over snow/ice) are 
assumed to have an ice fraction of 10/10. For regions 
identified as class 5 (broken sea ice) and class 7 (mixed 
surface types) the HHSC is used to estimate the frac- 
tional ice coverage from the visible data. Classes 1-3 
and 8-12 are assumed to contain exclusively land and 
water. The surface type beneath layered clouds is un- 
known in the algorithm and is treated as missing data 



Fig. 7. Automated cloud fraction analysis of the Arctic I July 1984 
imagery. Cloud fraction is given in tenths. 
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Fig. 8. Analyses of Northern Hemisphere (a) sea level pressure, 
and (b) 500 mb geopotential height at 12 UTC on 1 July 1984 
( ECMWF 1984). 


concentration throughout most of the region, with a 
mean absolute error of 0. 10. 

5. Algorithm characteristics 

This section examines the characteristics of the clas- 
sification and analysis algorithm. The assumptions used 
in the algorithm formulation are discussed, an analysis 
of the propagation of error is performed, and the hybrid 
histogram-spatial coherence method is compared to 
existing threshold methods. 

a. Discussion of assumptions 

The fundamental assumption underlying the entire 
approach is that an analysis ceil is homogeneous with 
respect to surface and cloud type. Examination ot the 
training imagery reveals that this assumption is gen- 



Fig. 9. Arctic surface weather and sky conditions observed on 1 
July 1984. The dashed line represents the 1/10 isopleth of sea ice 
concentration as analyzed by the NOAA-Navy Joint Ice Center tor 
3 July 1984. 


erally valid except in the presence of layered clouds or 
near coastlines. The former situation is ameliorated by 
the use of a category specifically meant for high clouds 
overlying low clouds. Class 7, mixed surface types, is 
intended to account for an inhomogeneous surface, 
but when clouds are present the classification is am- 
biguous. This is especially true when one surface is 
snow or ice and the other is land or water. In training, 
such a situation is classified as the appropriate cloud 
type over that surface type which appears to cover the 
most area in the cell. This problem is encountered often 
in the area of the Canadian Arctic Archipelago during 
summer. 



Fig. 10. Automated analysis of Arctic sea ice concentration tor 
1 July 1984. Values all given in tenths. 










May 1989 


ELIZABETH E. EBERT 


395 


Gouds, when present, are assumed to exist in a single 
layer. This is related to the previous assumption, and 
is one that is common to many, if not most, cloud 
analysis algorithms. The single layer assumption is fre- 
quently violated in the Arctic (Herman and Goody 
1976), and most severely near the core of cyclones. 
When two cloud layers are visible and well separated 
in the vertical domain, then class 18 (high cloud over 
low cloud) is the appropriate category. When one cloud 
layer completely covers a lower cloud layer, however, 
only the upper layer can be classified and analyzed. 
Situations of three or more layers of cloud also cannot 
be correctly analyzed by the algorithm. 

The algorithm assumes that the atmospheric long- 
wave contribution to upwelling radiances received by 
the satellite is negligible compared to the surface and 
cloud contributions. The water vapor content of the 
lower troposphere in the polar regions is typically quite 
low. The neglect of atmospheric water vapor absorption 
on the measurement of surface and cloud temperature 
results in an error of less than 0.5 K at 3.7 ^m and less 
than 2 K at 1 1 ^m (Duggin and Saunders 1984). Arctic 
haze, an aerosol attributed to anthropogenic sources, 
has its minimum concentration during the summer 
months (Rahn 1981). Over new snow, the presence of 
haze decreases the albedo measured from space, while 
over old snow and other darker surfaces, the albedo is 
slightly increased (Shine et al. 1984). These effects are 
not explicitly increased in the analysis algorithm. 

b. Error analysis 

In this section we examine the effects of classification 
errors. The magnitudes of the errors in cloud fraction. 


and surface and cloud albedo and brightness temper- 
ature resulting from cell m ^classification are largely 
dependent on the nature of the misclassification. Table 
7 illustrates the effect of several types of contrived clas- 
sification errors on the subsequent analysis. In these 
examples, incorrect classes were assigned to samples 
from the training set which were originally classified 
correctly by the algorithm. 

Certain classification errors may have little or no 
effect on the HHSC analysis if the confusion occurs 
between clouds with similar heights and underlying 
surface (type S). In this case the class-optimal channel 
does not change and the analyzed cloud fraction, clear 
and cloudy albedos and brightness temperatures are 
not significantly affected. 

When the misclassification involves a cloud of the 
wrong height (type H), the algorithm chooses an in- 
appropriate channel for the histogram analysis which 
may not show enough separation between the surface 
and cloud peaks. If the class characteristic value for 
the incorrect cloud type is used, as it was in the third 
example shown in Table 7, then the cloud fraction can 
be greatly under- or overestimated. 

A related error is the incorrect diagnosis of surface 
type underlying a cloud (type U). The analyzed surface 
albedo and brightness temperature are likely to be in 
error because the algorithm looks for a surface peak in 
the histogram within a reasonable distance of the (in- 
correct) class characteristic value. Finding none, it as- 
signs the albedo or brightness temperature to be that 
of the class characteristic value. Again, this can lead to 
large under- or overestimates in the fractional cloud 
cover. 


Table 7. Some examples of contrived classification errors and the resulting errors in the surface and cloud analysis. 

(T - true classification; F = false classification) 


Error type Class 

T = St over snow/ice 
F = Sc over snow/ice 
S 

T = Ci/Cc/Ac over land/water 
F - Cs/As over iand/water 

T = St over snow/ice 
F = Cs/As over snow/ice 

T = Ci/Cc/Ac over land/water 
F = Sc over land/water 

T = St over snow/ice 
F = St over land/water 

U 

T - Ci/Cc/Ac over land/water 
F = Ci/Cc/Ac over snow/ice 

T = St over snow/ice 
F = High cloud over low cloud 

T = High cloud over low cloud 
F = Cs over snow/ice 

T - St over snow/ice 
F - Unbroken sea ice 
N 

T - Cu over iand/water 
F = Land 


N 


T s 

ot c 

T c 

0.54 

0.56 

269 

0.62 

267 

0.54 

0.56 

269 

0.62 

267 

0.30 

0.10 

282 

0.42 

272 

0.30 

0.10 

282 

0.42 

272 

0.54 

0.56 

269 

0.62 

267 

0.03 

0.62 

269 

0.69 

251 

0.30 

0.10 

282 

0.42 

272 

0.55 

0.09 

282 

0.25 

275 

0.54 

0.56 

269 

0.62 

267 

0.98 

0.10 

271 

0.61 

268 

0.30 

0.10 

282 

0.42 

272 

0.09 

0.10 

282 

0.22 

251 

0.70 

0.61 

267 

0.64 

257 

0.88 

0.56 

267 

0.63 

253 

0.93 

0.10 

267 

0.46 

263 

0.59 

0.27 

265 

0.39 

258 

0.70 

0.61 

267 

0.64 

257 

0.00 

0.62 

254 



0.39 

0.10 

287 

0.23 

279 

0.00 

0.10 

286 
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The majority of misclassifications involved the con- 
fusion of “high cloud over low cloud” with single cloud 
types (type L). For clouds over snow and ice, errors 
associated with type L misclassifications result from 
the misinterpretation of the underlying surface as a 
lower cloud layer, thus overestimating the total cloud 
fraction, or from the misinterpretation of a lower cloud 
layer as snow or ice, in which case the cloud fraction 
is underestimated. If the region is south of 75° latitude, 
a misclassification of “high cloud over low cloud” may 
produce greatly erroneous estimates of surface albedo 
and brightness temperature if the class characteristic 
values are used (refer to the second type L example in 
Table 7). Over land and water, type L errors may be 
associated with multiple layers of the same cloud type. 
The total cloud fraction is not greatly affected in this 
case, but the partitioning of cloud fraction between 
low and high clouds will be in error. 

The most serious errors resulting from misclassifi- 
cation can occur when clouds in a cell are not detected 
at all (type N). Over land or water, the contrast between 
cloud and surface is usually great enough that only 
small amounts of cloudiness go undetected, but over 
snow and ice larger amounts of thin stratus or cirrus 
may be incorrectly analyzed as cloud-free sea ice or ice 
sheet. (This is problematic not only for satellite cloud 
observations, but also for surface observations of 
cloudiness during the polar night.) 

Of the 15% of the training samples which were mis- 
classified. 17% were of type S. 11% were type H, 2.o 
were type U, 52% were type L. and 18% were type N. 
For the contrived examples in Table 7, the errors in 
cloud fraction, surface and cloud albedo and brightness 
temperature associated with misclassifications were 
frequently quite large. These examples probably over- 
estimate the actual effects of misclassifications when 
they occur, however, because these scenes do not nec- 
essarily resemble the (incorrect) classes to which they 
were assigned. In section 3 it was shown that the sen- 
sitivity of the analysis to actual classification errors was 
not great, at least tor the analysis ot total cloud fraction. 

c. Comparison to threshold methods 

The hvbnd histogram-spatial coherence method of 
analyzing cloud fraction was chosen for two reasons. 
It takes advantage of the surface-cloud contrast in a 
class-optimal channel, and it accounts tor partially 
cloudy fields of view. This method is more complicated 
than the simpler threshold approach, and takes 25% 
more computing time per cell. It is worthwhile to com- 
pare the cloud fraction estimated using the HHSC to 
the threshold-estimated cloud fraction to evaluate 
whether the increased complexity is warranted. 

Three tvpes of threshold algorithms are tested. The 
first. (A), uses a VIS-IR threshold without classification. 
This is the most elementary type ot threshold method. 
A pixel must be both darker than the visible threshold 


and warmer than the infrared threshold to be labeled 
as surface; otherwise it is labeled as cloudy. This is the 
reverse of the usual logic, but yields better results at 
these latitudes than assuming a cell is clear unless 
proven cloudy. A step change in the threshold value is 
specified at 75° latitude to account for the surface tran- 
sition between an ocean and a sea ice regime. The sec- 
ond threshold algorithm, (B), classifies the cell first, 
then applies a VIS-IR threshold as before, with the val- 
ues of the thresholds now dependent on the class. In 
the third threshold algorithm, (C), the classification is 
followed by a threshold analysis of cloud fraction in 
the class-optimal channel only. In all three threshold 
algorithms, the values of the thresholds were computed 
from the class characteristic values and bounds (Ta- 

ble 3). L . 

Table 8 compares the results of the threshold anal- 
yses and the HHSC for the training set. The HHSC 
estimated the cloud fraction more accurately than the 
threshold methods in nine out of eleven cloud classes. 
The cirrostratus (classes 8 and 1 3) fraction was most 
accurately computed using (B), the VIS-IR threshold 
algorithm with classification. The HHSC was notice- 
ably more accurate than the threshold methods when 
analyzing low cloudiness, particularly over snow and 
ice. The greater success of the HHSC is due in part to 
its flexibility. It does not use a fixed threshold, but in- 
stead partitions the clear and cloudy pixels according 
to the statistical distribution of albedos or brightness 
temperatures in that particular scene. Classification 
improved the skill of the threshold algorithm in 
all cases. The one-channel class-optimal threshold 
method. (C), was more accurate than the VIS-IR 
threshold method. (B), for cirrus, stratus and strato- 
cumulus cloudiness. 

The improvement in cloud analysis achieved by first 
classifying a scene is exemplified by comparing the 
cloud fraction estimated by the HHSC and threshold 
(A) algorithm for the clear samples in Table 8. The 
threshold algorithm predicted cloud cover in a majority 
of the clear samples. Some of the clear samples actually 
do contain a very small amount of cloud (less than 
5%), but the contribution to the threshold RMS error 
is small compared to the gross errors incurred by mis- 
taking surface for cloud. Large overpredictions oc- 
curred for the snow on land, sea ice. and ice sheet 
classes which tended to be cold and bright in the im- 
agery. In contrast, the HHSC predicted cloud cover in 
less than 1% of the clear samples, the result of mis- 
classifications. The associated RMS errors are much 
smaller than those produced by the threshold algorithm 
without classification. 


6. Summary and discussion 

This study has demonstrated the detection and 
classification o! clouds in the polar regions trom multi- 
spectral satellite data using pattern recognition tech- 
niques. The algorithm is applied to 2.5° X 2.5° latitude- 
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Table 8. Comparison of the RMS errors for the hybrid histogram-spatial coherence method (HHSC) and threshold methods for the 
analysis of cloud fraction in the training set. Refer to Table 2 for category numbers. 






* 

RMS error 
Threshold methods 


Class 

Number of 
samples 

HHSC 

(A) VIS-IR without 
classification 

(B) VIS-IR with 
classification 

(C) Gass-optimal 
channel 

Surface 

1 

42 

0.00 

0.55 

0.00 

0.00 

2 

18 

0.08 

0.85 

0.25 

0.08 

3 

7 

0.00 

0.11 

0.00 

0.00 

4 

21 

0.16 

0.36 

0.15 

0.22 

5 

23 

0.16 

0.79 

0.03 

0.21 

6 

36 

0.09 

0.93 

0.15 

0.11 

7 

37 

0.00 

0.79 

0.00 

0.00 

Cloud 

8 

17 

0.34 

0.17 

0.11 

0.45 

9 

34 

0.25 

0.60 

0.42 

0.34 

10 

129 

0.16 

0.39 

0.24 

0.24 

11 

60 

0.15 

0.34 

0.35 

0.32 

12 

45 

0.11 

0.27 

0.25 

0.24 

13 

12 

0.24 

0.24 

0.23 

0.23 

14 

51 

0.35 

0.56 

0.49 

0.48 

15 

98 

0.26 

0.61 

0.59 

0.32 

16 

48 

0.18 

0.43 

0.44 

0.21 

17 

14 

0.06 

0.52 

0.35 

0.46 

18 

178 

0.23 

0.42 

0.39 

0.49 


longitude boxes of muitispectral AVHRR GAC-reso- 
lution data. Following the classification, the scene is 
analyzed to estimate cloud fraction, surface tempera- 
ture and albedo, and cloudy albedo and brightness 
temperature using a hybrid histogram-spatial coherence 
technique. When tested against manual estimates of 
cloud fraction, the mean absolute error was between 
0.05 and 0.26 for the training samples. The HHSC 
performed well when tested on sixty artificial scenes 
which had known values of fractional cloudiness and 
surface and cloud properties. 

One day of summertime satellite imagery from the 
Arctic was analyzed in detail. Large amounts of low 
and layered cloudiness were seen in the I July 1984 
analysis. The analyzed cloud types and fractions qual- 
itatively agree with manual interpretation of the sat- 
ellite imagery and with surface weather observations 
on the same day. Good correspondence between the 
algorithm's sea ice fraction estimates and the Navy- 
NOAA sea ice observations further validates the clas- 
sification and analysis scheme. 

The sensitivity of the classification and analysis to 
underlying assumptions and the propagation of error 
were explored in order to assess the abilities and lim- 
itations of the algorithm. Errors in classification can 
lead to large errors in the analyzed cloud fraction if the 
cloud height is incorrect, the incorrect surface type is 
diagnosed beneath the cloud, or if the cloud is not de- 
tected at all. Overall, the HHSC analyzes the fractional 
cloudiness more accurately than threshold methods 
using visible and infrared radiances or even one using 


the class-optimal channel. Comparisons between 
threshold methods with and without scene classification 
clearly show how the prior classification of surface and 
cloud type can greatly improve the accuracy of the 
cloud fraction analysis. 

This work demonstrates the ability of pattern rec- 
ognition techniques to successfully identify and analyze 
clouds in the polar regions, an area of the globe which 
confounds other cloud detection schemes. An algo- 
rithm such as this can be used to gather a data base of 
cloud observations in the polar regions, from which a 
cloud climatology can be constructed or numerical 
models can be validated. This is the goal of the Inter- 
national Satellite Cloud Climatology Project on a global 
scale. 

Several improvements may be made to the classifi- 
cation and analysis algorithm. This study was merely 
a demonstration of the use of pattern recognition tech- 
niques to aid in the analysis of cloud cover in a par- 
ticular region of the globe. Before such a scheme could 
be used operationally, more thorough training and 
testing would be required, including the use of a more 
extensive training set and a team of experts to train 
and evaluate the classification algorithm. The class 
characteristic albedos and brightness temperatures 
would be further refined in the process. 

The algorithm can be used to analyze regions which 
are smaller than the 2.5° X 2.5° chosen by the 1SCCP. 
Smaller samples would better satisfy the homogeneity 
and single cloud layer assumptions discussed in section 
5. Tests with the classification algorithm indicate that 
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its skill is not greatly reduced for training samples with 
dimensions of 16 X 16 pixels or even 8X8 pixels. The 
histogram analysis, however, requires at least 500-1000 
data points for accurate determination of the peak val- 
ues and widths. At 60° latitude, 500 GAC pixels would 
correspond to a sample with dimensions 1 .0° latitude 
by 1.25° longitude; at 80° latitude the sample dimen- 
sions would need to be at least 2° X 4°. 

Currently the higher resolution AVHRR data (1.1 
km field of view) is not saved over most of the polar 
regions, except by prior request. Use of this higher res- 
olution imagery would not only allow the minimum 
sample size to be reduced by a factor of 15, but the 
increased spatial detail in this imagery would almost 
certainly lead to more accurate classifications and cloud 
analyses using the pattern recognition and HHSC 
techniques. Tests are needed to determine the im- 
provement gained by using this higher resolution data, 
as well as the increase in computing resources that 
would be needed to perform the analysis. 

Temporal and spatial information could be utilized 
more effectively to improve both the classification and 
analysis. For example, the classifications of the nearest 
neighbors in both time and space might serve as use 
useful features for the classification of a cell, thus in- 
cluding some contextual information. The analyzed 
surface albedos and brightness temperatures of the 
nearest neighbors could also be used to filter out in- 
accurate analyses, at least over homogeneous regions 
such as open water, pack ice, or ice sheet. In cloud- 
free regions, cells classified as broken sea ice might be 
analyzed in more detail by subdividing the cell and 
employing a simple analysis algorithm (such as a 
threshold method) which does not require a large 
number of pixels. The higher resolution AVHRR data 
would be especially appropriate in this regard. 

Finally, additional data from the MSU 2 and HIRS 
instruments on the TOVS package could provide in- 
formation about the atmospheric temperature struc- 
ture, which would aid in cloud detection and height 
assignment. Microwave sea ice concentration mea- 
surements from the DMSP. SMMR. or SSM/I would 
be extremely useful tor specitying surtace characteris- 
tics. 

The next task, after analyzing the cloud type and 
amount, and surtace and cloud albedos and brightness 
temperatures, should be the estimation ot the cloud 
altitude and optical thickness and surtace radiative 
fluxes using radiative transfer models. The incorpo- 
ration of vertical temperature profiles would aid in the 
estimation of the atmospheric radiative fluxes. Class- 


: Kev to acronvms: DMSP. Defense Meteorological Satellite Pro- 
gram: HIRS. High Resolution Infrared Radiation Sounder. MSU. 
Microwave Sounding Unit: SMM/1. Special Sensor Microwave/Im- 
ager. SMMR. Scanning Multichannel Microwave Radiometer. TOVS. 
TIROS Operational Vertical Sounder. 


dependent microphysical models could be used to cre- 
ate a look-up table of cloud optical properties as a 
function of solar zenith angle, viewing angle, and cloud 
optical thickness. Before this can be done, however, 
much more needs to be learned about the microstruc- 
ture of clouds at high latitudes. Ideally this would in- 
clude both field measurements of cloud liquid and ice 
water content, drop size distribution, and crystal habit, 
as well as the development of improved methods to 
sense cloud microphysical properties remotely from 
satellites (e.g., Twomey and Seton 1980: Wu 1987; Platt 
etal. 1987). 
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Abstract 


A pattern recognition algorithm is demonstrated which classifies eighteen surface and cloud 
types in high latitude AVHRR imagery based on several spectral and textural features, then 
estimates the cloud properties (fractional coverage, albedo, and brightness temperature) using a 
hybrid histogram-spatial coherence technique. The summertime version of the algorithm uses 
both visible and infrared data (AVHRR channels 1-4), while the wintertime version uses only 
infrared data (AVHRR channels 3-5). 

Three days of low resolution AVHRR imagery from the Arctic and Antarctic during 
January and July 1984 were analyzed for cloud type and fractional coverage. The analysis 
showed significant amounts of high cloudiness in the Arctic during one day in winter. The 
Antarctic summer scene was characterized by heavy cloud cover in the southern ocean and 
relatively clear conditions in the continental interior. A large region of extremely low 
brightness temperatures in East Antarctica during winter suggests the presence of polar 

stratospheric cloud. 
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L Introduction and background 

The distribution of cloudiness in the polar regions is currently poorly known. This is 
largely due to the paucity of surface observing stations in the Arctic and Antarctic, but also due 
to the inherent difficulties in interpreting the much larger volume of satellite data collected in 
those regions. Certain characteristics of the polar regions make the analysis of cloud cover 
from satellite data more difficult than at lower latitudes, including low radiance contrasts 
between clouds and the underlying surface, the frequent inversion structure of the lower 
troposphere, and the lack of visible data during the winter. The key to improving our 
understanding of cloud cover, particularly in sparsely observed areas like the Arctic and 
Antarctic, lies in improving and increasing the use of satellite data for measuring standard 
quantities such as fractional coverage, cloud height and temperature, and if possible, cloud 
radiative and microphysical properties. At high latitudes, sensors aboard the NOAA series and 
Defense Meteorological Satellite Program (DMSP) polar orbiting satellites collect high 
resolution imagery which can be analyzed to estimate some of the above mentioned cloud 
parameters. 

Knowledge of polar cloud distribution is important for real-time applications such as 
weather analysis and forecasting, and is potentially valuable in initializing the moisture fields of 
numerical weather prediction models. Clouds control the hydrological and radiative processes 
linking the surface and atmosphere (Crane and Barcy 1984). Sensitivity tests conducted using 
general circulation models suggest that the wanning produced by an increase in atmospheric 
COt (the "greenhouse" effect) might be experienced most strongly in the polar regions (e.g. 
Hansen et al. 1984). Changes in the polar cloud distribution would likely accompany such a 
wanning trend. Regular monitoring of polar cloud cover might be able to detect this trend and 
its influence on the radiation balance via modifications in cloud conditions (WMO, 1988). 

In 1986 the International Satellite Cloud Climatology Project (ISCCP) launched a pilot 
study to investigate whether the information contained in the NOAA Advanced Very High 
Resolution Radiometer (AVHRR) channels could be effectively used in cloud algorithms to 
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resolve some of the cloud detection problems encountered in the polar regions. The AVHRR 
measures radiances in five spectral regions: two in the visible region (ch.l: 0.55-0.68 fim; 
ch.2: 0.725-1.1 fim), one in the 3.7 pm window region which is particularly useful for 
identifying low clouds over snow and ice during daylight (ch.3: 3.55-3.93 (im), and two 
thermal infrared window channels (ch.4: 10.5-1 1.3 pm; ch.5: 11.5-12.5 pm). A pilot data set 
comprising fourteen days of Global Area Coverage (-4 km resolution) AVHRR imagery from 
January and July 1984 was analyzed by various research groups using a variety of schemes, 
and the preliminary results were presented in a workshop report (WMO, 1987). 

Trained observers of satellite imagery identify surface and cloud types at high latitudes by 
recognizing such features as cracks in the sea ice, the smooth homogeneity of ice sheets, partial 
obscuring of surface details by thin clouds, cloud shadows and illumination of cloud sides, and 
the bumpy or wispy appearance, as well as the usually cooler brightness temperatures of 
certain cloud types (McGuffie et al. 1988). Automated pattern recognition uses the same 
textural and spectral (radiance) features to objectively classify the components of a scene. The 
textural features in particular are extremely useful in distinguishing clouds from sea ice and ice 
sheet in a low contrast scene. Ebert (1987a) used ten spectral and textural features in low- 
resolution visible, near-infrared, and infrared AVHRR imagery to classify eighteen surface and 
cloud types in Arctic summertime data using the maximum likelihood method. Welch et al. 
(1989) found that textural measures in a single visible channel of very high resolution 
LANDSAT imagery were sufficient lo classify four surface and cloud types in tire Arctic, again 
using the maximum likelihood method. Key (1989) used several spectral and textural 
features to identify cloud patterns in Arctic AVHRR imagery. Yamanouchi et al. 
(1987) used brightness temperature differences in AVHRR channels 3, 4, and 5 to 
detect, but not classify, clouds over Antarctica. 

Other schemes have been developed to identify’ clouds over snow and ice from spectral data 
alone (Bolle 1985, Molders 1987, Gesell 1989, Key et al. 1989a.b). Spectral methods classify 
each pixel individually, producing a classification image with the same resolution as the 
original imagery, while methods employing textural measures operate on clusters of pixels and 
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generally produce a lower resolution classification image. 

The majority of polar cloud schemes have been developed and tested on daytime data, 
where cloud identification relies largely on reflectance contrasts in the visible and particularly in 
the near- infrared portions of the solar spectrum. During the dark half of the year, only thermal 
radiances are available for use in cloud analysis. Of the methods mentioned above, only 
Yamanouchi’s algorithm is not affected by the loss of visible data. 

This paper presents some results from a pattern recognition analysis of a portion of the 
ISCCP polar cloud pilot data set. The pattern recognition algorithm of Ebeit (1987a) is adapted 
for wintertime data using AVHRR channels 3, 4, and 5. The cloud types and fractional 
coverage in 2.5° latitude x 2.5° longitude regions (corresponding to the ISCCP grid) are 
analyzed in three days of imagery from Arctic winter, Antarctic summer, and Antarctic winter. 
The results are discussed in the context of concurrent synoptic weather observations. 


2. Method 

The classification and analysis scheme will be only summarized here; a full description is 
given in Ebert (1987a,b, 1989). The algorithm consists of two parts: (1) the classification of 
surface and cloud type in each 2.5° x 2.5° scene using pattern recognition, and (2) the analysis 
of cloud fraction, clear and cloudy visible albedos and infrared brightness temperatures using a 
hybrid histogram-spatial coherence technique. 

The purpose of the classification step is to determine whether a scene contains clouds, and 
if so, to identify the general physical nature of the clouds. It does this by computing several 
spectral and textural features from the calibrated multispectral AVHRR data, then using these 
feature values in the maximum likelihood method to classify the scene into one of eighteen 
surface and cloud categories. A different combination of AVHRR channels and classes is 
necessary for each season and hemisphere. The Arctic summer and Antarctic summer versions 
of the algorithm use AVHRR visible and infrared channels 1, 2, 3, and 4, while the Arctic 
winter and Antarctic winter versions use the infrared channels 3, 4, and 5. For the Arctic 
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winter algorithm the categories are: 


1. Land without snow 

2. Land with snow 

3. Open water and thin ice 

4. Unbroken sea ice 

5. Broken sea ice 

6. Ice sheet 

7. Mixed surface types 

8. Cs/As over land/water 

9. Ci/Ci/Ac over land/water 


10. St over land/water 

1 1. Sc over land/water 

12. Cu over land/water 

13. Cs/As over snow/ice 

14. Ci/Cc/Ac over snow/ice 

15. St over snow/ice 

16. Sc over snow/ice 

17. Inversion St over snow/ice 

18. High cloud over low cloud 


(The standard abbreviations for cloud type are from the International Cloud Atlas , WMO 1956) 
In the two summertime algorithms class 3 was restricted to "open water” and class 17 was 
assigned to "Cu over snow/ice" to represent cumulus and other isolated or small-scale low 
cloudiness. No large areas of snow-free land exist south of 60 S, so class 1 was reassigned 
to "ice haze" in the Antarctic summer algorithm and "polar stratospheric cloud" in the Antarctic 
winter algorithm. 

An iterative training procedure was used to teach each classification algorithm the 
characteristics of the surface and cloud categories, based on knowledge of the correct 
classification of a training set of image data. Approximately 1000 samples were included in 
each training set. The algorithm's classification ability was evaluated using the skill score, S, 
defined by 

_ # correct classifications - # expected by chance ^ ^qq 
# total classifications - # expected by chance 

The classification skills achieved on the training data are shown in Table 1 . The greatest skill 
was achieved for the two summertime versions, reflecting the usefulness of visible imagery. 
The Arctic winter algorithm had the least skill, due mainly to the difficulty in detecting low 
clouds over snow and ice from thermal data alone. All of the versions had some problems 
separating layered cloudiness from single layer clouds. Tests suggest that the classification 
skill on independent data is a few percent less than that achieved on the training data (Ebert, 
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1987a). 

Overall the presence of cloud in a scene was correctly detected 97% of the time during 
summer and 89% of the time during winter. Figure 1 shows a confusion matrix of predicted 
versus actual presence of cloudiness for each of the four versions of the pattern recognition 
algorithm. The majority of clear/cloudy errors resulted from cloudy regions being mistaken for 
clear surface, especially in winter. 

The classification of a scene specifies the optimal AVHRR channel for separating clear and 
cloudy pixels in the hybrid histogram-spatial coherence technique. For example, channel 4 
separates clear and cloudy pixels best for "Ci/Cc/Ac over land/water," while channel 1 is 
optimal for "Cu over land/water" and channel 3 is optimal for "St over snow/ice" during 
summer. The surface and cloud classifications are also used to specify initial guesses for the 
clear and cloudy radiances. The incorporation of the classification information is an extremely 
important step in the cloud analysis process. Ebert (1989) showed that the prior classification 
of a scene can greatly improve the determination of cloud fraction in the polar regions where 
the surface-cloud contrast may be extremely small in the visible and/or infrared imagery. 

The hybrid histogram-spatial coherence technique estimates the total cloud fraction, and the 
clear and cloudy channel 1 albedos and channel 4 brightness temperatures. Briefly, a 
histogram analysis is first performed in the class-optimal channel to estimate the mean clear 
(|i s ) and cloudy (|i c ) albedos or brightness temperatures. Next, the spatial coherence method 
(Coakley and Bretherton, 1982) isolates completely clear and completely cloudy arrays of 2x2 
pixels, which are tagged and given weights of w=0 or w=l, respectively, with the weights 
representing pixel cloud fraction. Partly cloudy arrays are given weights w = (x-p. s )/(p. c -p. s ), 
subject to the constraint that 0<w£l. The average of all the weights is an estimate for the total 
cloud fraction. The partly cloudy pixels are then removed, and the completely clear and 
completely cloudy pixels are used to generate "clean” histograms in the remaining spectral 
channels, from which the clear and cloudy albedos and brightness temperatures can be 


estimated. 
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When verified against manual nephanalyses of the training imagery, the overall root-mean- 
square (RMS) error in the estimated cloud fraction was 0.18 for the Arctic summer samples 
and 0.24 for the Arctic winter samples. The classification step ensured that most of the clear 
scenes were correctly identified. The summertime algorithm analyzed the low cloud cover well 
with a RMS error of 0.18, but had more difficulty analyzing the high cloud cover, with a RMS 
error of 0.31. The Arctic wintertime algorithm performed better on the high cloudiness (RMS 
error of 0.23), but underestimated the low cloud cover (RMS error of 0.35) because some low 
level clouds were not detected. The Antarctic summer and winter algorithms were not tested 
against manual analyses. However, Fig. 1 suggests that the accuracy of the cloud fraction 
estimated by the Antarctic algorithms do not differ greatly from those of their Arctic 
counterparts. 

3. Analyses 

Each day's satellite data comprises fourteen images covering most of the region between 
60° and 90° latitude. The images were analyzed individually, then the analyses combined and 
averaged to form a hemispheric mosaic. The region immediately surrounding the pole was not 
analyzed because those pixels, found at the edges of the imagery, were distorted due to the 
earth's curvature. In this section the analyses are compared to manual interpretation of the 
imagery and to the cloud cover observed at the surface. 

3.1 6 January 1984 Arctic 

Figure 2a shows a mosaic of thermal (channel 4) imagery for the Arctic on 6 January 1984. 
The European Centre for Medium Range Weather Forecasting (ECMWF) analysis of sea level 
pressure at 12 UTC, as well as surface synoptic observations of cloud cover from the 
European and Canadian Arctic (Free University of Berlin 1984; NOAA 3-hourly surface 
weather analyses) at 7 UTC and 15 UTC, respectively, are plotted in Fig. 2b. In using the 
surface observations to validate the automated analysis, it should be noted that these 
observations are point measurements made at a specific time, and may not represent the daily 
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mean conditions over a larger area which is portrayed by the automated analysis. 

The imagery shows an extended region of cumulus cloudiness over the Norwegian Sea, 
with thick low clouds over northern Sweden and Norway. Finland and western Russia are 
partially clear under the influence of a ridge. Layered and high clouds are found in Siberia 
between 60°N and 65°N. The Barents Sea between Svalbard and Novaya Zemlya contains 
clouds at all levels, over open water in the south and sea ice in the north. The Canadian 
Archipelago and Baffin Bay are mostly clear, with broken and unbroken sea ice and some 
regions of warm stratus. The textural variability indicates the presence of some cloud cover 
over Greenland, but it is difficult to tell whether there is cold high cloudiness over a wanner 
surface, or warm stratus over a colder surface. The Arctic Basin appears to be mostly clear 
with some warmer inversion stratus in the eastern region. 

The classification and cloud fraction analysis of the image are shown in Figs. 2c and 2d. 
The algorithm accurately classified the cumulus clouds in the Norwegian Sea, the clear 
conditions over the Arctic Basin and the Canadian Archipelago, and the high cloudiness in 
eastern Greenland. The high clouds near 90°E were also well classified. It failed to correctly 
diagnose the thick low cloud cover over northern Scandinavia, and the presence of sea ice east 
of Novaya Zemlya. It is difficult to verify the clouds analyzed over central Greenland. 

3.2 6 January 1984 Antarctic 

A mosaic of visible (channel 1) imagery from 6 January 1984 over Antarctica is shown in 
Fig. 3a. Only one satellite analysis was made for some regions north of 70°S due to the poor 
illumination, as the satellite descending orbit passed over during the early morning hours. The 
homogeneous appearance of the ice sheet suggests clear skies over much of the interior. The 
ocean surrounding the continent is mostly ice-free and quite cloudy, with particularly bright 
cloud cover over the western Weddell Sea near 50°W. Figure 3b shows the surface analysis 
from January 6. The map shows three major cyclones situated at 60° latitude, with the 
strongest of them located near 10°W longitude. The clouds to the east of the Antarctic 
Peninsula may be related to convergence at the mountain barrier (Schwerdtfeger 1984). An 
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anticyclone is situated over East Antarctica, with few clouds there. 

The classification of the image and analyzed cloud fraction are presented in Figs. 3c and 
3d. Stratus and stratocumulus were the most frequent cloud types analyzed over the oceans. 
Layered cloudiness is located primarily in the eastern sides of the cyclonic systems. The 
clouds over the ice sheet are mainly of liquid phase, and are found throughout the lower 
elevation regions of Antarctica. The high plateau is mostly clear, with areas of warm, tenuous 
ice haze near 20°E and 100°E. A strong gradient in cloudiness is observed at the coast near 
90°E, the possible result of cold drainage flow from the ice sheet acquiring moisture over the 
warm water and condensing to form clouds. The analyzed cloud types do not agree very well 
with the small number of available surface reports, which indicate alto and cirriform cloudiness 
where the algorithm analyzes low cloudiness. 

3.3 1 July 1984 Antarctic 

The infrared imagery for the 1 July 1984 Antarctic is combined into a mosaic in Fig. 4a. 
The accompanying surface synoptic analysis is given in Fig. 4b. East Antarctica appears 
homogeneous but the measured 11 pm brightness temperatures (not shown) range from 179- 
200 K, as much as 19° colder than the surface temperature observed at Vostok. The quality of 
the AVHRR calibration at temperatures below 200 K is not well known, but the errors in 
brightness temperature are probably not more than 1-2°C (see Appendix A). These very cold 
temperatures seem to signal a broad region of polar stratospheric cloud. As in the summertime, 
several cyclonic disturbances surround the continent, with low and layered cloudiness near 
their centers and eastern flanks. Cloud-free sea ice is apparent east of the Antarctic Peninsula 
and again between 60°E and 90°E. High cloudiness covers the coast at the eastern side of the 
Weddell Sea, and thin cirrus covers much of West Antarctica. 

The classification and cloud fraction analysis of the imagery are presented in Figs. 4c and 
4d. Unbroken sea ice was analyzed in and east of the Weddell Sea and in the region near 
60°E. East Antarctica was mostly clear, except for the two large inland lobes of polar 
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stratospheric cloud near 40°E and 100°E. Small amounts of cirrus and warm inversion stratus 
were analyzed over West Antarctica. Inversion stratus was analyzed mainly around the margin 
of the continent. There is fair agreement between the analyzed cloud cover and the surface 
observations. Unfortunately it is not possible to confirm the presence of the analyzed polar 
stratospheric cloud using limb scanning measurements, which depend upon the extinction of 
sunlight (McCormick et al. 1982). 

4. Conclusions 

The analyzed cloud fields show many of the features which have been recorded m 
climatologies of the Arctic and Antarctic (e.g. Vowinckel and Orvig 1970; Schwerdtfeger 
1984), such as the wintertime cumulus clouds over the Norwegian Sea and high cloudiness 
over Siberia, the continual presence of multiple levels of clouds accompanying cyclones in the 
southern ocean, and the predominantly clear interior of the Antarctic continent during summer. 
No climatology exists for polar stratospheric clouds, whose existence has been known for 
barely a decade (McCormick et al. 1982). These clouds are found in both hemispheres during 
winter when the stratospheric temperature falls below about 195 K, and may play a catalytic 
role in the springtime destruction of Arctic and Antarctic ozone (Hamill et al. 1986). 

Automated pattern recognition has been shown to be an effective tool for analyzing cloud 
cover from satellite data in the polar regions, an area where many conventional techniques 
encounter difficulty. The pattern recognition algorithm demonstrated here classified 2.5 x 
2.5° scenes with a skill of 71-86%, depending on the season and hemisphere, and estimated 
the cloud fraction to an accuracy of 0.18-0.24 as measured by the RMS error. While these 
results are encouraging, they fall far short of the accuracy of ±0.03 for 30-day mean cloud 
fraction required by the ISCCP (Schiffer and Rossow 1983). Some improvements could be 
made by decreasing the dimensions of the analysis region from 2.5° x 2.5° to a smaller size 
which would better satisfy the homogeneity assumptions implicit in the algorithm, and by 
increasing the number and quality of the training samples to improve the classification skill. 



10 


Better results would almost certainly be obtained using the higher resolution ( 1 . 1 km) AVHRR 
data in which the finer scale features, particularly of low level clouds, are more evident. 
Microwave and infrared radiometry and high spectral resolution interferometry (Smith et al. 
1988) give additional information on the surface conditions and atmospheric structure which 
would aid in cloud detection and height assignment. With the improvement of technology to 
remotely sense and analyze polar clouds, we can continue to increase our knowledge of their 
distribution and improve our understanding of the role they play in the global climate. 
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Appendix A. Low temperature calibration of AVHRR thermal sensors 

The radiance, R, measured in an AVHRR thermal channel is estimated as a linear function 
of the raw counts, X: 

R = S X + 1. 

S and I are the slope and intercept values valid for each scan line, determined by a 2-point fit to 
a set of in-flight calibrations against an internal black body target and outer space. 

The radiance is converted to brightness temperature by inverting the Planck function in the 
appropriate spectral region. Hie dependence of the Planck function itself on temperature is 
accounted for by varying the central wavelength, X, of each band over three temperature 
ranges. For channel 4 of the NOAA-7 AVHRR these central wavelengths are X=10.797 pm 
for T= 180-225 K, X=10.790 pm for T=225-275 K, and *=10.785 pm for T=275-320 K 

(Kidwell, 1984). 

The 2-point linear calibration does not account for the nonlinearity of the sensor response 
over the observed range of radiances. This can result in temperature errors of up to 3°C 
(Brown et al. 1985). A recent NOAA publication (NOAA 1986) lists values of the nonlinear 
correction term which should be added to the linearly calibrated brightness temperature to 
obtain a more accurate value of brightness temperature. These corrections, plotted in Fig. 5, 
were determined from an analysis of the original thermal vacuum data which did not include 
measurements at temperatures below 205 K. However, the trend clearly favors negative 
corrections at colder temperatures. A possible source of error is contamination of the scene 
radiance by internal emission from the radiometer housing (Brown, personal communication). 

The AVHRR data in this paper were calibrated and analyzed prior to the publication of the 
NOAA memorandum using only the linear calibration outlined by Kidwell (1984). Figure 5 
suggests that these temperatures may be 1-2°C warmer than the true scene brightness 

temperatures. 

A post-launch calibration of thermal channel radiances against known values of scene 
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temperature (for example, hourly surface temperature measurements made at several Antarctic 
stations) might be used to test the linearity of the radiometer response at temperature below 205 
K. In the future all pre-flight calibrations of radiometers intended for operation on polar 
orbiting satellites should include some very low temperature measurements, as these are not 
uncommon at high latitudes during winter. 
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Classification 
skill score 


Problem areas 


Arctic summer 

85.8% 

Layered clouds, small amts, of Cu over land 

Arctic winter 

71.9% 

Layered clouds, low clouds over snow/ice, 



especially inversion St 

Antarctic summer 

85.1% 

Layered clouds, some Ci over ice sheet 

Antarctic winter 

78.4% 

Layered clouds, low clouds over snow/ice, 



especially inversion St 


Table 1. Classification skill scores achieved on four training sets. 


Figure captions 


Fig. 1. Confusion matrices of clear and cloudy classifications for the (a) Arctic summer, (b) 
Arctic winter, (c) Antarctic summer, and (d) Antarctic winter training samples. 

Fig. 2a. AVHRR channel 4 image mosaic for the 6 January 1984 Arctic. (The routine for 
producing the image mosaic uses a maximum brightness criterion, and thus the mosaic appears 
colder and cloudier than the individual images due to cloud motion between the fourteen 
satellite overpasses.) 

Fig. 2b. ECMWF sea level pressure analysis at 12 UTC and observed cloud conditions on 6 
January 1984. 

Fig. 2c. Automated surface and cloud classification for the 6 January 1984 Arctic. (Key: 
l=Land without snow, 2=Land with snow, 3= Open water and thin ice, 4=Unbroken sea ice, 
5=Broken sea ice, 6=Ice sheet, 7=Mixed surface types, 8=Cs/As over land/water, 9=Ci/Cc/Ac 
over land/water, 10=St over land/water, ll=Sc over land/water, 12=Cu over land/water, 
13=Cs/As over snow/ice, 14=Ci/Cc/Ac over snow/ice, 15=St over snow/ice, 16=Sc over 
snow/ice, 17=Inversion St over snow/ice, 18=High cloud over low cloud) 


Fig. 2d. Automated cloud fraction analysis for the 6 January 1984 Arctic. Cloud fraction is 
expressed in units of tenths. 

Fig. 3a. AVHRR channel 1 image mosaic for the 6 January 1984 Antarctic. 

Fig. 3b. ECMWF sea level pressure analysis at 12 UTC and observed cloud conditions on 6 
January 1984. 



Fig. 3c. Automated surface and cloud classification for the 6 January 1984 Antarctic. (Key: 
l=Ice haze , 2=Land with snow, 3=Open water, 4=Unbroken sea ice, 5=Broken sea ice, 6=Ice 
sheet, 7=Mixed surface types, 8=Cs/As over land/water, 9=Ci/Cc/Ac over land/water, 10=St 
over land/water, ll=Sc over land/water, 12=Cu over land/water, 13=Cs/As over snow/ice, 
14=Ci/Cc/Ac over snow/ice, 15=St over snow/ice, 16=Sc over snow/ice, 17=Cu over 
snow/ice, 18=High cloud over low cloud) 

Fig. 3d. Automated cloud fraction analysis for the 6 January 1984 Antarctic. Cloud fraction is 
expressed in units of tenths. 

Fig. 4a. AVHRR channel 4 image mosaic for the 1 July 1984 Antarctic. 

Fig.4b. ECMWF sea level pressure analysis at 12 UTC and observed cloud conditions on 1 
July 1984. 

Fig. 4c. Automated surface and cloud classification for the 1 July 1984 Antarctic. (Key: 
l=Polar stratospheric cloud, 2=Land with snow, 3=Open water and thin ice, 4=Unbroken sea 
ice, 5=Broken sea ice, 6=Ice sheet, 7=Mixed surface types, 8=Cs/As over land/water, 
9=Ci/Cc/Ac over land/water, 10=St over land/water, ll=Sc over land/water, 12=Cu over 
land/water, 13=Cs/As over snow/ice, 14=Ci/Cc/Ac over snow/ice, 15=St over snow/ice, 
16=Sc over snow/ice, 17=Inversion St over snow/ice, 18=High cloud over low cloud) 


Fig. 4d. Automated cloud fraction analysis for the 1 July 1984 Antarctic. Cloud fraction is 
expressed in units of tenths. 

Figure 5. Nonlinear temperature correction as a function of the channel 4 brightness 
temperature obtained using the 2-point linear calibration for the NOAA-7 AVHRR (from 


NOAA 1988). 
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